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LARGE DEVIATION PRINCIPLE FOR INTERACTING BROWNIAN 

MOTIONS 


INSUK SEO 


Abstract. We prove the Large Deviation Principle for the empirical process in a system of in¬ 
teracting Brownian motions with singular interactions in the nonequilibrium dynamic. Such a 
phenomenon has been proven only for two lattice systems: the symmetric simple exclusion process 
and zero-range process. Therefore, we have achieved the third result in this context and moreover 
the first result for the diffusion-type interacting particle system. 
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1. Introduction and Outline 

1.1. Introduction. The large scale behavior of tagged particles is a primary concern in interacting 
particle systems. The first breakthrough was accomplished by Kipnis and Varadhan m, whose 
seminal paper introduced a general invariance principle for additive functionals of reversible Markov 
processes. Furthermore, these authors derived the equilibrium central limit theorem (CLT) for the 
tagged particle in the symmetric simple exclusion process (SSEP) as an application of the general 
theory. This equilibrium result has been extended to various models, e.g., interacting diffusions 
[21| . the mean-zero asymmetric simple exclusion process (ASEP) by Varadhan |35) and the ASEP 
for d > 3 by Sethuraman, Varadhan and Yau |29| . 

The study of tagged particles in the nonequilibrium is a field with many untapped possibilities. 
Recently, several researchers have developed nonequilibrium CLTs for tagged particles especially for 
ID interacting particle systems. Jara and Landim |14) proved such a result for a ID nearest neighbor 
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SSEP where the equilibrium CLT had been established by Rost and Vares m- Sethuraman and 
Varadhan |28) investigated the corresponding large deviation principle (LDP). The nonequilibrium 
CLT for the ID systems has also been proven for the nearest neighbor ASEP by Goncalves [9], 
simple exclusion process with long jumps by Jara SSI, zero-range process (ZRP) by Jara, Landim 
and Sethuraman |15| and locally interacting Brownian motions by Grigorescu |10| . 

Another approach is to study the empirical process Rn = which can be regarded as 

the averaged tagged particle. In the 1990s, Quastel, Rezakhanlou and Varadhan found a systematic 
approach to study the limit theory and large deviation theory for the empirical process, which 
resulted in a series of published work |23l|2ll[26l[33]. Despite the robustness of their methodology, 
the LDP for this context is only known for two models: the SSEP in the case where d > 2 and 
the ZRP. The limited applicability of their general method arises from the lack of intermediate 
large deviation theory for the empirical density of colors and such a theory was only available for 
the SSEP |24| and ZRP jT]- We refer to the survey paper by Varadhan m for a comprehensive 
discussion of this research. 

The main purpose of the current work is to explore the third result for this context. Our inter¬ 
acting particle system is the locally interacting Brownian motions on the one-dimensional torus T 
that were introduced by Grigorescu mm- The main result of |10| is the nonequilibrium CLT and 
asymptotic independence of two tagged particles for an interacting Brownian system and together 
these fulfill the law of large numbers (LLN) for the empirical process. However, the approach that 
was followed in |10| does not rely on the empirical density of colors, but on the hydrodynamical 
analysis of the local time and hence the LDP of the kind proposed in |24) was unavailable. In this 
study, we have established the LDP by analyzing the empirical density of colors. 


1.2. Interacting Diffusion with Local Interaction. We start by introducing the interacting 
Brownian motions on T with local interactions^ In this model, we assume that N Brownian 
particle^ • • • , Xjy(-) are moving on T with partial reflections. To state this succinctly, two 

particles reflect each other when they collide, but sometimes they change their labels. We can 
measure the collision time between two particles xf{-) and {■) up to time t by the local time and 
then the switching of labels occurs as a Poisson process with the constant rate \N{X > 0) along this 
canonical local time. We now provide a rigorous definition of this informally described interacting 
particle system and introduce the notion of empirical density of colors. 


1.2.1. Definition of Particle System. Our particle system x^ft) = [x^ ft), ft), ■ ■ ■ , xfff{t)) can 
be regarded as a diffusion process on T'^. We can define this diffusion in three equivalent ways: 
generator, martingale problem formulation, and Dirichlet form. All of them are of course useful 
during our excursion. 

^For the detailed definition of model, we refer |101 lll| . 

^We also able to assume that the number of particle is ajv such that limjv_>oo ^ = p for some p > 0. 
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We first introduce an A^-manifold 

Gn = {x = (ti, X 2 , ■ ■ ■ , xn) £ : Xi 7 ^ xj for all i / j} 

and consider the boundary OGn = £>i<j{xi = xj}. Note that each face {x : Xi = xj} consists of 
two sides. We will follow the convention that Fij is the side at which Xj approaches Xi from the 
clockwise direction so that with the usual orientation of T, Xj = Xj + 0 on Fij whereas xj = Xj — 0 on 
Fji- We only consider piecewise smooth functions on Gat and that are smooth up to the boundary 
such that we can define 

fij{x) = /(••• , Xi-i, Xi + 0, Xi+i, • • • , Xj-i, Xj - 0, Xj+i, ■■■) 

Dijfix) = (Vi - Vj)/(- • • , Xi_i Xi + 0, Xi+i, • • • , Xj-i, Xj - 0, Xj+i, ■■■) 

for X G Fij for all i, j. Let us denote this class of functions by G{Gn). We are now in a position to 
define the process x^{t). 


(1) Generator: The generator for the process is = 5^/ s-iicl the domain P(jSf7v) consists 
of functions / G C{Gn) that satisfies the boundary condition iifjf{x) = 0 on Fij for each 
i ^ j where 

^ijf{x) = Dijf{x) - \N{fij{x) - fji{x)). (1.1) 

This infinitesimal generator {.^n, defines the process. 

(2) Martingale problem formulation: We define the process on by the measure Ptv on 

(^([O, T], T'^) for some fixed final time T. Under this measure, we have N{N — 1) local 
times and the filtration {^t : t G [0, T]}, such that for any f £ C (Gn), 

Mf{t) =/(!'"(()) - /(^^''(O)) 



is a martingale with respect to{^t:0<t<T}. The martingale Mf{t) can be represented 
in another way such that 

^ rt f‘t 

= '^kf{x^{s))dj3k{s)+ '^ {fji - fij){x^{s))dMjJ(s) (1.3) 

i^jJo 

where /3k{t), 1 < k < N is a family of independent Brownian motions and 

M!J{t) = j!^{t)-XNAf^{t) 

where J^j , l<iy^j<N is a family of pairwise orthogonal Poisson jump processes with 
rates XNAfj{t); hence, (t) as well as (t)^ — XNAfj{t) are martingales for each i j. 
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(3) Dirichlet Form: For each / E C{Gm), the Dirichlet form is given by 

^N{f) = 11 |V/(x)| ^dx + ^Yl [ - fji{x)fdSij{x) (1.4) 

^ JGn ^ JFij 

where dSij{x) is the Lebesgue measure on Fij normalized to have the total measure to be 
l|^ It should be noted that the first part of (ll.4p corresponds to the Brownian movement of 
particles, whereas the second part takes into account the Poisson jump type of interaction. 
Now, we have the process x^{t) on T'^ for t E [0, T] which can be regarded as a system of diffusion 
processes with local interactions. 

Remark 1.1. Since is a covering space of T^, we can lift any continuous trajectory in to the 
one in in a unique fashion. Therefore, we implicitly regard x^{t) as a process on sometimes. 
This does not cause any technical issues, because we usually work with the density field of the form 
Tf f i^)) with a periodic function / on M. For the detailed explanation, see Section 1.4 of 

uni. 

To understand the large scale behavior of our interacting particle system, we start by considering 
the empirical density given by 

1 ^ 

= (‘-5) 

i=l 

for 0 < t < T which induces a measure Qtv on (^([O, T], For the level of the empirical 

density, our process is equivalent to the non-interacting case, because our interactions essentially 
involve the switching of labels which does not affect (II.Sp . The limit theory and large deviation 
theory of the empirical density in the non-interacting system is well known. To state such results 
in a concrete form, we need the following assumption. 

Assumption 1. The initial empirical density satisfies the LLN in the sense that 

/r'^(O) ^ p^{dx) weakly in ^{T) for some non-negative measure p^{dx) on T with a total mass of 
1. Moreover, also satisfies the LDP with the rate function ) o-nd scale N. 


Then the LLN and related LDP of the empirical density can be formulated as follows. 


Theorem 1.2. Under Assumption 1, converges weakly to the Dirac mass on the trajectory 

{p{t, x)dx : 0 < t < T} where p{t, x) is the solution of the heat equation dtp = with initial 
condition p^{dx). Moreover satisfies the LDP with the rate function 


1 /■’ 

■)) + 2 J 


da - 2 ^^ 


dt 


( 1 . 6 ) 


-1,7 


and scale N where the H ^ norm is defined in the standard way. 


^Note that the usual normalization for the Lebesgue measure on this diagonal face is U2 
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Proof. See mM for the LLN and LDP, respectively. 


□ 


Remark 1.3. The LLN and LDP for the interacting system of diffusive particles have been developed 
at the level of the empirical density for several models. The first results were published in the classic 
papers mm for the LLN and LDP for weakly interacting Brownian motions on M, respectively. In 
this model, the interaction comes into play through the drift coefficient of the form 6(xj(f), 
Recently, |22l [2] developed corresponding theories for case in which particles interact through their 
ranks and hence the diffusion coefficients also depend on In both models, the limiting 

dynamics are governed by McKean-Vlasov type of equations. Another type of result was presented 
in |34) for the two-body model 

dxf {t) = —N ^ V'{N{xf {t) — xf{t))dt + dj3i{t) (1.7) 

on T, where V{-) is a compactly supported, even and smooth potential. In this model, the limiting 
particle density is given by the unique solution of dtp{t, x) = ^[P(p(f, x))]xx where P{-) is the 
pressure functional depending on V. Even if the author did not explicitly establish the LDP, 
estimates therein are sufficient to establish the LDP through standard methodology of [1]. It was 
already pointed out in |10| that our local interaction model can be regarded as a limit of the 
two-body interaction model. More precisely, if we consider the sequence of the potential {I4}e>o 
satisfying 

lim / exp{2I4(a:) — l}dx = --do 

in the sense of distribution, then the corresponding diffusion given by (Il.Tjl with I4(') converges 
to our locally interacting diffusion P^r with parameter A. We refer the readers to |10| for details. We 
also remark here that the totally asymmetric counterpart of our model has been studied in UM- 


1.2.2. Empirical Density of Colors. We now introduce the notion of the empirical density of colors, 
which is an intermediate object toward the empirical process. 

Let I 2 , • • • ) dm} be a (non-random) partition of{l, 2, such that satisfies 


tN\ 


lim 


= Pc >0 


( 1 . 8 ) 


N^oo N 

for each 1 < c < m where pc is the average density of color c. Then, {xf^{■) : i E 1^} denotes the 
set of particles of color c and empirical density of color c is defined as 


1 


V(‘) = E 


Ad 


(1.9) 


ieip 


for t E [0, T]. Finally, the empirical density of colors is defined by 

= id-lit), d2it), ••• , R 




( 1 . 10 ) 
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which induces a probability measure Qjsr on (^([O, T], In contrast to the uncolored em¬ 

pirical density (|1.5I) . we have to take the interaction into account, because the switching of labels 
between particles of different colors affects in a complex manner. Similar to the uncolored 

empirical density, we need an assumption on (0) to obtain the limit theory and large deviation 
theory. 


Assumption 2. The initial empirical density of colors satisfies the LLN in the sense that 

fP{dx) = (/5?(dx), pl{dx), • • • , pl,{dx))^ (1.11) 

weakly in where the non-negative measure p^{dx) has a total mass of Pc for each c. More¬ 
over, also satisfies the LDP with the rate function and scale N. 

Note that Assumption 2 implies Assumption 1 with p^ = 

1.3. Main Results. Based on the precise description of the model in the previous section, we now 
summarize our main results. 


1.3.1. Large Deviation Theory for Empirical Density of Colors. The hydrodynamic limit theory for 
{Qiv}?^=i is implied by the results of |10| . It is well known [261 [3T] that if the scaling limit of the 
tagged particle is the diffusion with the generator TCp, which may depend on the limiting particle 
density p, and any two tagged particles are asymptotically independent, then the limiting particle 
density pc of each color is the unique weak solution of the PDE dtPc = .^p Pc with the initial 
condition p^{dx) for each color c. 

The scaling limit as well as the asymptotic independence of the tagged particle of our system 
has been studied m- The scaling limit turned out to be the diffusion with the time-dependent 
generator 


— _ 

^ 2{X + p{t,x)) 


^ _ (2A x))Vp{t, x) ^ 


( 1 . 12 ) 


2{X + p{t, x))2 

where p{t, x) is the solution of the heat equation as in Theorem 11.21 Consequently, under Assump¬ 
tion 2 , Pc is the solution of the parabolic equation 


dpc 

dt 



X 

X-\- p 


V Pc + 


Vp 
X-\- p 


(1.13) 


with the initial condition p^{dx) for each color c. We are also able to reorganize these equations in 
the form of a matrix with the notation p = (pi, p 2 , ■ ■ ■ , Pm)^ as 

^ = ^V-[D(p)Vp] (1.14) 

with initial condition p^{dx). The m x m diffusion matrix D{p) is explicitly given by 


6ij X-\- Pi 
X-\- p 


D{p)ij 


(1.15) 
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In other words, Qtv converges weakly to a Dirac mass on the unique solution of (|1.14p . We explain 
the details with another proof of this result in Section 4. 

We can decompose the diffusion matrix by D[p) = A[p)x{p) where A{p) and x{p) defined by 




^ij^Pj 4“ PiPj 
X + p 


X{p) = diag 


i- - 

I ’ ’ 

\Pl P2 



(1.16) 


Here, x{p) is the Hessian of the entropy functional h{p) = P* iogp* and A{p) = D{p)x{p)~^ ■ 
Note that this matrix A[p) is symmetric, which is not a coincidence; rather, it is a so-called Onsager 
reciprocity (see [7] for details). 

Under Assumption 2 and an additional technical assumption related to the uncolored initial 
profile (Assumption 3 in Section 4.1), the LDP for can be established with the rate 

function L™;o^(t.) for vr. G (^([O, T], ^(T™)). More precisely, L™;o^(vr.) < oo only if Tit is absolutely 
continuous with respect to the Lebesgue measure for each t, and if we write such vr. as p{-, x)dx 
then 


ImloriK-, x)dx) = Cit(p(0, x)dx) + lXyn{p{-, x)dx) 
where /j™jj(-) and L^„(-) explain the large deviation rates of the initial configuration and the dynamic 
evolution of the system, respectively. The initial rate function is just a part of Assumption 

2. The dynamical rate function L^„(-) is our primary concern and is given by 


I^nipi-, x)dx) = 


-f 

2 Jo 


dt 2 


dt. 


(1.17) 


[Dip)Vp] 

The rigorous meaning of this expression is carefully explained in Section 4.2. This large deviation 
result is the main contribution of the current work and is explained throughout Sections 2, 3 and 
4. In Section 2, we establish some super-exponential estimates which essentially mollify the local 
times into the local densities. Section 3 provides the exponential tightness of {Qa}^= 1 ) which 
compactify the upper bound problem. Relying on these preliminary results, we compute the exact 
upper and lower bounds in Section 4. 


1.3.2. Large Deviation Theory of Empirical Process. The empirical process is defined by ^ ^xP{- 

which induces a measure P/v on ^i(C'([0, T], T)). As we mentioned earlier, P/v converges weakly 
to the delta mass on a diffusion P G .J^i(C(lO, T], T)) of which the generator is given by jj/p defined 
in ()1.12p . Our main concern is the LDP corresponding to this result and developed in Section 5. 

We now explain the main result. For a probability measure Q on (^([O, T], T), we can explain 
the LDP via the rate function ^{Q), which is finite only if Q has marginal densities at any time 
t G [0, T]. Let us denote this marginal density by q{t, x), then ^{Q) < oo only if q is weakly 
differentiable in x and satisfies 


/ (?(0, x) log g(0, x)dx < oo and 
Jt 


f IVgP 

/t Q 


■dxdt < oo. 
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For such Q, we consider a class of function on [0, T] x T given by 

= |6(f, t) : dtq = ^Aq - V [bq], IL b^qdxdt < oo 

where the sense of PDF is weak. Then we can find the unique diffusion process with the generator 
^q,b = + bV and starting measure q{0, x) for each b G ^q. Note that also has the marginal 

density q{t, x). Then we can prove that the relative entropy H{b) := H[Q\P^] is either finite for all 
b G or identically infinite. We set ^{Q) = oo for the latter case. For the former case, we can 
find a 6q G which minimizes H{-) on SSq, and then the rate function is given by 

rT 

' 

IQ JT 

The last two terms measure large deviation rates from P to P^'^ and to Q, respectively. 

The final important remark is that this LDP result is a direct consequence of that of the empirical 
density of colors. The LDP rate for the finite dimensional projection of the empirical process can be 
understood as the one for the empirical density of colors. By doing so, we can obtain the full LDP 
by using Dawson-Gartner’s projective limit theory. This profound relationship has been revealed by 
Quastel, Rezakhanlou and Varadhan in |24) and, therefore, our work not only verifies the robustness 
of their methodology, but also the universality of their large deviation result. This is because our 
rate function (ll.lSp is quite similar to that of the SSEP, whereas the model dynamics are seemingly 
unrelated. We will explain the robustness and universality in Section 5. 


J^{Q) = /mit(g(0, •)) + ^ / [ x)q{t, x)dxdt + H 

^ Jo Jt 


Q 


P^Q 


(1.18) 


2. Super-exponential Estimates 


In the LDP theory in the context of interacting particle systems, the core step is to establish the 
replacement lemma, named after Guo, Papanicolau and Varadhan’s seminal work |12| . which is a 
super-exponential type of estimate. We introduce an appropriate form of the replacement lemma 
and related concepts in Section 12.11 We prove this by first providing some preliminary estimates in 
Section 12.21 following which the proof of the main replacement lemma will be given in Section [ 


2.1. Replacement Lemma. In general, the empirical density is studied via its corresponding 
density field. In the context of the hydrodynamic limit theory as well as the large deviation theory 
for the empirical density, the main difficulty is the replacement of the current, which appeared at 
the computation of density field as a result of interactions, by gradients. Due to the nature of our 
dynamic, the current is highly related to the local time as we can see from (11.21) . Consequently, 
the replacement lemma should replace the local times by the appropriate gradient type, which we 
achieved by first introducing the averaged local times and local densities. 

The basic local times in our process are 


A?{t) = lim [‘ = /‘ - Pmdt 
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for i ^ j where 5+ is a delta-type distribution on [0, 1] such that f{x)6^{x) = ^/(O). We can 
understand as the amount of the collision time between two particles x^{■) and xf{-) up to 

the time t at which the particle xf{-) approaches to x^{') from the clockwise direction. Of course, 
each local time A^{t) is quite noisy and impossible to estimate alone. Fortunately, these noises can 
be controlled by taking average among them. Two such examples are 


•^'"(*) = V2E[A?(*) + A'i(‘)] 

i¥=3 

(2.1) 

T(‘) = L2 Kw + A'lw] 

(2.2) 


for each i. Note that Af {t) is the average collision time of the particle xf{-) against all the other 
particles up to time t and A^{t) is the total average of local times. The behavior of these averaged 
local times has been studied extensively in m and is also important for our work. However, as our 
focus is on the density field of the empirical density of colors the main object to be estimated 

is 

<=(*) = V E [A?W + PA 

which measures the average collision time of the particle xf{-) against the particles of color c. 

Now, we define the notion of local densities. For x = (xi, X 2 , • • • , xat) G T'^, the local density 
function of color c around the particle Xi is defined by 

Here, the function is the usual indicator function on T and we henceforth simply denote 

this function by Xe(')- The essence of the replacement lemma for our model is the replacement of 
the integral with respect to the average local time of the form 

dAf^M = {dAfjit) + dAfiit)) 

3&Ic 

by the usual integral of the form p^fl{x^{t))dt. Formally, this replacement can be stated as the 
following theorem. 


Theorem 2.1 (Replacement Lemma). For e, <5 > 0, 0 < < ^2 E T o-nd two colors ci, C 2 , let 

t 2 ',e, 6) be the event defined by 




* 64 ^ -^*1 






LARGE DEVIATION PRINCIPLE FOR INTERACTING BROWNIAN MOTIONS 


10 


Then, we have 

limsuplimsup — logP'^ t 2 ]e, 5)1 = —oo. (2-4) 

N^oo N 

The main object of the remaining part of current section is to prove this theorem. 

2.2. Preliminary Estimates. 


2.2.1. Green’s Formula for Gm- We frequently use Green’s formula on the A^-manifold G^- As 
a calculation of this nature is not conventional, we briefly explain our philosophy in this short 
subsection. 

Let us denote the zth standard unit vector by e*. Then, we can apply Green’s formula for 
f £ C{Gn) and vector field V(x) = with V) G C{Gn) for each i, so that 


[ (V/(x), V(x)) dx 

JGm 

= — f /(x) (V • V) (x)(ix + f /(x) (V(x), n(x)) (iS'(x) 

Jgn JOGn 

= - [ fix) (V • V) {x)dx + ^ / fij{x) (v(x), Ci - Bj) dSij{x). 
JGn TT JFii 


(2.5) 


Remark 2.2. Recall from (HID that dSij{x) is the Lebesgue measure on Fij normalized to have a 
total measure of 1. The unit normal vector of the boundary Fij is — ^j), but is eliminated 

from (12.51) because of this renormalization. 


In particular, some special forms of the vector helds provide us useful results. We summarize 
such results by the following lemma. 

Lemma 2.3. For f G C{Gn), h G ^^([0, 1]) and Vi{x) = Ylj-j^i h{xj — Xi), the vector field 

N 

V{x) = ^Vi{x)ei. 

satisfies 


2 = 1 


[ (V/(x), V(x))dx 

JGn 


/ fix) (V • V) {x)dx 
Jgn 

h(l) - h{0) 


yZ / + fjiix)] dSijix). 

In addition, suppose that Ui{x), • • • , C/Ar(x) G C(T'^) satisfy Uiix) = Uj{x) whenever Xj = Xj, for 
all i ^ j. Then, the vector field 

N 


W(x) = ^17,(x)R,(x)e, 


2 = 1 
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satisfies 


[ (V/(x), W(x)) dx 

Jgn 


/ fix) (V • W) {x)dx 
Jgm 

/i(l) -/i(0) 


y] / Ufix) [fijix) + fjiix)] dSijix). 


Proof. For the first part, it is enough to check boundary terms. Note that 

(Vix), ei - Cj) = hixj - Xi) - hixi - Xj) + [hixk - Xi) - hixk - xfi] 

k:k^i,j 

and hixk—Xi) — hixk—Xj) = 0 on Fij for k i, j. Moreover, hixj—xi) = /i(0) and hixi — Xj) = /i(l) 
on Fij and hence 


V /" {Vfix),Vix))dx 
i=i ■’Gn 

= - [ fix) (V • V) ix)dx - (h(l) - h(0)) V [ fijix)dSiji. 
JGn JFii 


'Gn i^j 

by (12.51) . Obviously, we can symmetrize the last term as 


X] [ fiji^)dSijix) = [ ifijix) + fjiix)) dSijix) 




2 Jp.. 


and we are done. The proof of the second part is similar. 


□ 


2.2.2. Estimates Based on Dirichlet Form. The proof of Theorem 12.II heavily relies on the Dirichlet 
form. In particular, we frequently use D^iVy) for some / > 0 and we denote this by V^if). By 

(m, 


= i L ^ ? L [v^ - 


dSijix). 


( 2 . 6 ) 


'Gn ^ V-/ ijLj 

In addition, let be the class of non-negative functions / € CiGip) which also satishes fj fix)dx = 

1 . 


Lemma 2.4. For any f G we have 


1 

N 


i=l •'Gn 

Proof. This bound can be proven by 
N 


N 

X] / |Vi/(x)|dx< 

. _ 1 J Gtkf 


SVNif) 
N ' 


(2.7) 




|V/(x)| dx < 


y/N 


iv/(^)r 

fix) 


dx / fix)dx 


□ 
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Lemma 2.5. For any f £ we have 


1 

F 2 


f ifijix) + fjiix)) dSijix) <2 + 


( 2 . 8 ) 


Proof. Define a function h{x) = x — ^ on [0, 1], then Vi{x) = ~ Xi) satisfies ViVi{x) = 

— {N — 1) and |Vi(x)| < {N — l)/2. Hence, by Lemma [2^ 

\Y ! dSijix) 


i¥^j 


JFi 


N{N - 1) / f{x)dx + Y [ 


< N{N-1) + 


N - 1 


TV 

V/ \VJix)\dx. 


We can complete the proof by Lemma 12.41 


□ 


Let us define M^^i{x) = Ylj jj^iXeixj — Xi) which counts the number of particles around Xj. The 
following series of lemmas provides estimates related to Me^j(-). We also remark here that we shall 
write C for a constant and as usual different occurrences of C may denote different constants. 


Lemma 2.6. For any f £ m o-nd e > 0, we have 


[ f{x)Y^d,4x)dx<C + 

JGn j=i \ / 


Ar 2 

Proof. We take an auxiliary function on [0, 1] such that 

for 0 < X < e 


(2.9) 


h^{x) = < 


X — € 


0 fore<x<l — e 

X — (1 — e) for 1 — e < X < 1 

Note that Vi(x) = Ylj-j^i^eixj — x,) satisfies |Vi(x)| < (N — l)e and ViV)(x) = —Me^j(x). Thus, 
by Lemma 12.31 

TV 

/ f{x)Y^Gii^)dx 

= Y [ Nifi^)]^ii^)dx + eY [ ifijix) + fjiix)) dSij{x) 

i=l 


i^j •’Fj 


a:) + fjiix)) dSij{x) 


<iN-l)eY[ |Vi/(x)|dx + e^ / {fij{ 

,=i Jgj, 

and therefore we can complete the proof by applying Lemma 12.41 and 12.51 


□ 
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Lemma 2.7. For any f £ ^]\f and e > 0, we have 

1 r ^ 

lv2 y] |Vi/(x)| < C 

Proof. By Cauchy-Schwarz’s inequality, 

1 




1 

€2 . 


( 2 . 10 ) 


^ [ y]|Vi/(x)|M,,i(x)(ix <^(f dx j f{x)J2Mliix)dx 

JGn^—i \JGn J\^) JGn j^—i J 

1 

N \ 2 


< 


1 8VNif) 1 


( JV 'N^ja„ 

since M^^i{x) < N. Thus, (12.101) is direct from Lemma 
Lemma 2.8. For any f G CLnd 0 < e < j, we have 


„ N \ 

/ f{x)'^M^^i{x)dx 

Jgn , 


□ 


1 r ^ 


ffs I ■/ zJ ^ ^ 

JGn i=i 






( 2 . 11 ) 


Proof Since M^^{x) = Ylk,j^iXe{xk - Xi)xe{xj - Xj), we need to bound Xe{x)xe{y) by a more 
tractable object. To this end, we define a function ke{-, •) G C'(T^) for e < |. Firstly, along the 
y — X = c with 0 < c < 2e, kf:{x, y) is defined by 


I line 


ke{x, y) = 


if —e<x<e — c 
jE§{^ — e — c — x) lie — c<x<l — e — c 


2e- 


0 


ifl — e — c<x<l — 


For y — X = c, —2e < c < 0, 


-i y — a. — ^ ^ u, we set k^{x, y) := k^{y, x). Finally, k^{x, y) = 0 for all the other x, y. 

It is easy to see that G C'(T^) and moreover k^ satisfies 

( 2 . 12 ) 

(2.13) 


0 < ke{x, y) < 2 ex 2 e{x - y) 

Xe{x)Xe{y) < ‘2ex2e{x - y) + Vxke{x, y) + Vyk^ix, y). 

In particular, (I2.13P enables us to bound M‘^^{x) such a way that 

Mli{x) < M^^i{x) + y] [2eX2e{Xp-Xq)-ViK{Xp-Xi,Xq-Xi)\ 

p,q:p,q^i,p^q 




N 

i(x) + e y M 2 e,i{x) - ViK^^i{x) 
i=i 


K^^i{x)= y keixp - Xi, Xq - Xi). 

p,q:p,q^i,p^q 


where 
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Therefore, we have 
N 


y] / f{x)Mli{x)dx<'^[ f{x)M^^i{x)dx + eN f f{x)'^M 2 e,i{x)dx 
^ Jgn Jgn Jgn 

N 

- / f{x)S/iK^^i{x)dx. (2.14) 

i=i Jgn 

We can bound the first two terms of the RHS by Lemma 12.61 For the last term, we can apply 
Green’s formula with the vector field K(x) = so that 


N N 

y^ / f{x)ViKi^^{x)dx = V' / Vif{x)Ki^^{x)dx. 
Jgj, Jgn 


(2.15) 


Note that boundary terms are disappeared since is continuous on T^. Moreover, by (|2.12p . 


N 


Ki^e{x) < 2e X2e{Xp - Xq) < e'^M2e,l{a 


(2.16) 


p¥=q 


1=1 


Consequently, we can bound the last term of (|2.14p by (12.151) and (|2.16p . such that 


yy f Vif{x)Ki^^{x)dx < ( f 

Gn \ J G 


iv/(x)r 


■dx 


N 

/ f{x)'^Kl^{x)d. 

JGm 


'Gn j=i 
N 


< 


Gn 

IGeN^VNif) [ f{x)^K,4x)d3 
Jgn 

r ^ 

VM / f{x) Ney^M2,,i{: 

Jgn 


< 16eN‘^ 


1=1 


dx 


where we used the trivial bound K^^i{x) < 2eN‘^ at the second inequality. Now, we can complete 
the proof by applying Lemma 12.61 □ 


Lemma 2.9. For any f ^ o-nd 0 < e < j, we have 




Me,i{x){fij + fji){x)dSij{x) < C 
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Proof. Let us define two auxiliary functions Ug, on [0, 1] by = Ji Xeiy)dy and tte(x) = 

fi v^{y)dy. They enjoy the following properties: 

u'^{x) = Ve{x) and n'(x) = Xeix) (2-17) 

e 2 e 

— Xf (a;) < Ue{x) < —Xe{x) and \Ve{x)\ < -Xe{x) (2.18) 

1 X^( 0 ) = Ue{l) = y and ^^(O) = -^,Ve{l) = (2.19) 

Let US denote V'e^j(x) = ~ Xi) and Ue^i{x) = — Xi). Then the vector field 

W(x) = U^^i{x)Ve^i{x)ei satisfies conditions of the second part of Lemma 12.31 and therefore, 

Pi = P2 + Ps (2.20) 


where 


-Pi = ^ f Ue,i{x) ifijix) + fjiix)) dSij{x) 


i¥=j 


N 


P 2 = - V' / f{x)Vi [Ue,i{x)Ve,i{x)] dx 
i=l JGn 
N 

P 3 = -V'/ U^,i{x)Ve,i{x)\/if{x)dx. 
i=l JGm 

We first estimate Pi. Notice that ^^(l) — ^^(O) = e by (12.191) and Ut^i{x) > j(x) by (12.181) . 
These together give ns 

Pi > X] / ,*(®) + fjiix)) dSij{x). (2.21) 


Now, let US consider P 2 . By (|2.17p . 

-Vi [P,,i(x)W,i(x)] = P,,i(x)M,,i(x) + V,^i(x), 

which is bounded above by |e^MP(x) due to (12.181) . Thus, we can bound P 2 as 


P 2 < CN^ 




\ N J 


7 e 


( 2 . 22 ) 


by Lemma 12.81 

We now bound P 3 . Since Ue^i{x)V^^i{x) < ^MP(x) by (12.181) 

N 


IP3I < 


E/ 1^‘A 

1 Gat 


if{x)\ M^ i{x)dx < CN'^ 




62 


i=l '’Gn 

by Lemma [2.71 Now, (I2.20p . (I2.2ip . (I2.22p and (|2.23p implies the desired bound 


(2.23) 

□ 
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The next and last preliminary estimate controls the discontinuity of / along the boundary, joint 
with Mj^e(x). 

Lemma 2.10. For any f G 0,^(1 0 < e < we have 


L “ fjiix)\Mi^^{x)dSij{x) < C 






\ N J 


1 1 
64 + —= 

Vn 


Proof. By Cauchy-Schwarz’s inequality, 


L \fijix) - fjiix)\Mi,eix)dSij{x) 


< 


1 


{\/“ \/ dSij{x) 

[ (\Jfijix) + \Jfjiix)) ^le{x)dSij{x) 


J F 

.A 21 


< 


1 /4 Pjv(/) 


^ / (/d(^) + fii{x)) Ml^{x)dSij{x) 




A^2 I A N 

and the proof is completed by Lemma iT 
2.2.3. Main Estimate. 

Proposition 2.11. For any f ^ and 0 < e < we have 

1 


Ar2 


E 


' Gn 


_ 1 

^ 2 e 2 


fij{x)dSij{x) + fji{x)dSji{x) 


< C 




i 1 

64 + —= 

^/n 


□ 


(2.24) 


Proof. Note first that the presence of absolute values in the summation prevents us from applying 
Green’s formula in the form of Lemma [2.31 Instead, we define a function on [0, 1] as 

X 1 


9eix) = < 


2e 2 
0 


for 0 < X < 6 
for 6 < X < 1 — 6 


^ + 1 for 1 - 6 < X < 1. 


and then consider a vector held V(x) = geixi — Xj)ej consisting of only one direction. Green’s 
formula for this vector held is 


V. 


'Gjv 


f{x)ge{Xi - Xj)dx = f f{x)^XeiXi - Xj)dx + Ji + J 2 
JGn 


(2.25) 
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because g'^{x) = ^x>l{x) where Ji and J 2 are the boundary terms to be explained below. The 
boundary terms for this particular Green’s formula are 


E 

p^q L 


/pq(3^) (V(x), Cp Gq) dSpqi^x) 


(2.26) 


for which the summands are non-zero only if p or gr is j. Now, let Ji be the sum of two summands 
in (I2.26|l with (p, q) = (i, j) or (j, i) and J 2 be the sum of all the others. Namely, 

Ji= fji{x)ge{xi-Xj)dSji{x)- fij{x)ge{xi - Xj)dSij{x) 

J Fji J Fij 




fjk{x)ge{xi - Xj)dSjk{x) - / fkj{x)ge{xi-Xj)dSkj{x) 


Note that g^{xi — Xj) = ^ on Fij and — ^ on Fji and hence 




1 


1 


Ji = -o fji{x)dSji{x) - - fij{x)dSij{x). 


I F ■ 

JTni 


fij{x)dSij{x) + fji{x)dSji{x) 


(2.27) 


For J 2 , we know that ge{xi — Xj) has same value on Fjk and Fkj and therefore 


E 




{fjk{x) - fkj{x))ge{xi - Xj)dSjk{x) 


(2.28) 


By combining (I2.27p . (|2.28p with (I2.25P we can bound the LHS of (I2.24p by J 3 -|- J 4 where 




[ \'^jf{x)\\ge{Xi- Xj)\dx 


= / \f3k{x) - fkj{x)\\ge{Xi- Xj)\dSjk{x) 

A-LA lr,.U—LA A " 


i^j k:k^i,j 


Since \ge{-)\ < \xe{-), we have 


N 


J 3 < / |Vi/(x)| M,j{x)dx 


J4 E 


2 iV 2 

1 


— 1 JGn 


i=l 


i^j jk 

The proof is completed by Lemmas 12.71 and 12.101 


Y [ \fijix) - fij{x)\ M^j{x)dSij{x) 

A —L A Fi U 


(2.29) 

(2.30) 

□ 


2.3. Proof of Replacement Lemma. In this subsection, we provide the proof of Theorem 12.11 
based on Proposition 12.11] and the classic technique developed by Donsker and Varadhan [1]. Their 
method is only available for the process Pat which is sufficiently close to the equilibrium process 
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in the sense that 


log 


(JFa 


L°°(T) 


= 0{N). Unfortunately, this condition does not hold not only 


for our model, but also for the general interacting particle system of diffusion type. For example, 
if we start deterministically, our process Pjv is even orthogonal to P^ starting from the invariant 
measure dx. We solve this issue by using a symmetrization procedure. However, this procedure is 
only possible for the time slot [rj, T] for some r] > 0. Thus, we have to establish the replacement 
lemma on the interval [0, g] in an independent manner. Let us examine this procedure more closely, 
by dividing the Theorem 12.11 into the following two propositions. 


Proposition 2.12. For any rj, 6 > 0, r] < ti < t 2 < T and two colors ci ^ C 2 , 

limsuplimsuplogPAT t 2 ]e, 5)1 = —oo. 

Af^.oo ^ 

Proposition 2.13. For any 5 > 0 and two colors ci ^ C 2 , 

lim sup lim sup lim sup 
77^0 e^O N^oo 



(2.31) 


(2.32) 


In this paper, estimates in the form of Propositions 12.12] and [27131 are referred to as by normal time 
regime and small time regime, respectively. These respective regimes require different approaches. 
Dichotomies of this nature frequently occur in our work. 


2.3.1. Normal Time Regime. The key procedure on which to base the proof of Proposition 12.121 
is symmetrization. Our interacting particle system is defined as a probability measure Pat on 
(^([O, r], T'^) with the initial profile f^{dx); then, the density profile at time t > 0 denoted by 
/Ar(t, x) satisfies the forward equation 

x) = ^Afjv(t, x) + x)S+(xj - Xi) (2.33) 

i¥=j 

where ii^/ is as defined in (jl.ljl . The process with the initial density dx is the equilibrium process, 
which we denote by P^. 

Now, we define some intermediate processes. Let be the set of all permutations of [N] and let 
a{x) = (Xo-(i), T(t( 2 )j • • • j 3:o-(Ar)) for cr G TPa^ and x G Then we can consider a process starting 
from 

fNidx) = ^ fN{da{x)) (2.34) 

o-e*Pjv 

with the same interacting mechanism. We denote this process by Pat. Finally, we define another 
initial profile 




(2.35) 
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where C is the set of all permutations with • • • , (cf. Section 1.1.2) as their 

invariant sets. Then, let be the process with the initial profile with the same type 

of interactions. 

Lemma 2.14. Let Eat he an event on (^([O, T], T^) which only depends on sub-path {t(s) : rj < s < 
T}. Furthermore, if the event Eat is invariant under permutations in the sense that {t(-) G Eat} = 
{(t(x(-)) G Eat} for all a G ^Pa^; then we have 

N 

(2.36) 


Piv [Ejv] < P^ [Ejv] 


for some universal constant C. Furthermore, if the event Eat is only invariant under the permu¬ 
tations among the same color in the sense that {x(-) G Eat} = {(t(x(-)) G Eat} for all a G (Ln, 
then 


Pat [Eat] < 


Cm 


N 




(2.37) 


where m is the number of eolors. 

Proof. First of all, the marginal density of the process Pat is 

hit, a;) = ^ X] 

o-eTiv 

and hence, we can deduce from (I2.33h that hit, x) is the solution of the heat equation dtfN = ^A/at 
with initial condition (12.341) . Therefore, we have bound of the form 


\\hit^ OIL ^ 


_c 

Vi 


N 


<-ivr 


(2.38) 


for at > r] with a (universal) constant cB Note that if Eat is invariant under all permutations then 
Ptv [Eat] = Pv [Eat]- Moreover, since Eat only depends on the path after time rj, we have 


Pa [Eat] < 


Vn 


N 


Pa"[Ea] 


by (I2.38P and therefore we can derive (|2.36p . 

For (j2.37l) , note first that the marginal density prohle of P“^°^ at time t is 

1 


fV^Vt, x) = 




hit, crix)). 


o-eCjv 


Since |eiAf| = A^i!A^ 2 ! • • • where = jL |, we can obtain 

fV°Vt, x) _ ^ fNit, Vx)) 


< 


N\ 


fNit,x) fNit, crix)) NilN 2 l---Nm 


< m 


A 


^It is easy to check since fN{t, x) = f% * Pn{x) where p%{x), the heat kernel on T, is given by 
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and thus OHoo — t > r] from (I2.38p . Therefore, we can derive (I2.37P in a 

similar way. □ 


If r/ < < t 2 < T, then the event t 2 ;e, S) satisfies the conditions of second part of 

previous lemma and thus 


Fat (f)] < 


Cm 

Vv 


N 


F?,"[C"^’"^(ti,t2;e, <5)]. 


Consequently, we can reduce Proposition 12.121 into the following equilibrium estimate. 


(2.39) 


Proposition 2.15. For any r], 5 > 0, r] < ti < t 2 < T and two colors ci ^ C 2 , 

limsuplimsup ^ logP^ t 2 ]e, <5)] = -oo. 

£^>0 N^oo ty 


Proof. By Chebyshev’s inequality, 

^logP^ t2\e, 5)] 


< -ad + — log exp 


{•sif 


{t))dt - K>2(t2) - 


for any a > 0 where E^ is the expectation with respect to P^. Note that 

E^^expio^ f it))dt - {Af^^^{t2) - 

[ *64^ 

< E^expja^ei ^ pff (x^ it))dt - {A^C2it2) - A^^^ih)) 


Ci=±lVi 


and let us investigate each summand of the last line. By Feynman-Kac’s formula, 

pt2 
fti 

< exp{(t 2 - ti)XN,e,a} 

where A 7 v,e,a is the largest eigenvalue of the operator 


E^^exp I a ^ e* ^ pif{x^{t))dt - (A^ih) - AYiti)) 




(2.40) 


(2.41) 


(2.42) 


(2.43) 
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on The variational formula for is 


sup I ^ I 


[ /(^ 

Jgn 


XejXj - Xj) 

2e 


dx 


1 


fij{x)dSij{x) + / fji{x)dSji{x) 


-VnU) 


By Proposition 12.111 the expression inside sup (and thus \N,e,a) can be bounded by 


N 




£4 + 


y/N 


VnU) 


N 




< C’N 

where C, C are proper constants. Therefore, by (I2.43p . (|2.42p is bounded by 

2'^exp{C'A^(ae4 + + OAr(l))(t2 - ti)} 


(2.44) 


Finally, by (I2.4ip and (I2.44h . 

^logP^ [C^’‘'2(ti, t 2 ;e, d)] < -a6 + log2 + C{ae^ +o®e^)(t 2 - h) + OAr(l) 
and therefore 

limsup limsup ^ log[C^’'^^(ti, t 2 ;e, (5)] < -a(i + log2. 

e^O TV^-oo ^ 


Since a > 0 is arbitrary, we are done. 


□ 


Remark 2.16. Previous argument combining Chebyshev’s inequality, Feynman-Kac’s formula and 
the variational formula for the maximal eigenvalue, in the context of interacting particle system, has 
been originally introduced by |3] and also explained thoroughly in Chapter 10 of m- This method 
will be used frequently and implicitly in the remaining part of the current article. 


2.3.2. Small Time Regime. We now prove Proposition 12.131 By Chebyshev’s inequality. 


N 


log Pat 


N 


E 


- Afm,) 


> 


< —a5 + — log Eat exp 

< —a5 + — log Eat exp 


E 


p[f {x^ {t))dt - 


N 


N 

pe,i{x^ {t))dt + aN (rj) 


or any a > 0 where Eat is the expectation with respect to Pat. Accordingly, we only need to establish 
the following estimates. 
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Proposition 2.17. For any a > 0, 


lim sup lim sup — log Eat exp { aN ( 17 )} < 0 

Tj^O N—>-oo 


lim sup lim sup lim sup — log E n exp 
r/^O e^O N—>-oo 



< 0 


(2.45) 

(2.46) 


Proof. This proposition is a special case in which the labels of particles (and thus the interaction) 
play no role. Therefore, for this proposition we can temporarily assume that x)^(t), X 2 {t), ■ ■ ■ , x^{t) 
move by the way of independent Brownian motions. For x E T'^, let us dehne G]\f{x) = ^ 

Xj) where g{x) = is a continuous function on T. Then, by Tanaka’s formula, 

N( \\ , Nfnw , ar]{N - 1) 2a 


Gn{x^^v))-Gn{x^^0)) + 


N 


-V r 


Y1 9 \xf{s)-xf{s)) 


N 


dxi{s) 




i¥=j 


and thus by Girsanov’s theorem, 

Eat exp {2aNA^{r]) - [GAr(x^(r/)) - GAr(x^(0))] - An{v)} 
aT]{N — 1) 


where 


= exp ■ 


A7v(i/) = 


2 iV 2 


E/ E/(-f( 

*=1 “ \pi^i 


s) — xf (s)) I ds < 


afigN 


since |fl'^(T)| < By (I2.47h and (I2.48P we can obtain 


EArexp{2aA^A^(?7) - [GNix^ir])) - GAr(x^(0))] } 
< exp 


[ arjN a^rjN 


Moreover, by the mean value theorem, 

Eat exp { Gat (t^ (??))-Gat (x^( 0 ))} < Eat 


N 


exp 


iv) - xf {0)\ 


< < 1 + a 


2=1 

TT 


exp 


2 

at] 


N 


(2.47) 


(2.48) 


(2.49) 


(2.50) 
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since we have assumed that xf {t), 1 <i < N are independent Brownian motions. By (|2.49p . (|2.50p 
and Cauchy-Schwarz’s inequality, 


Eat exp {a( 77 )} < 



ar]N 5a^riN^ 

^ 16 j 


(2.51) 


and we proved (|2.45p . 

For (j2.46p . we define Pe{x) = ^Ue{x) where Ue is the function defined in Lemma 12.91 and thus 
p”{x) = ^Xe{x)- Define Hn{x) = ^ P^i^i ~ ^j) ^-PPly Tanaka’s formula such that 


HM{x^'^{p))-HM{x^m-2a 


N 


Pe,i {x^{t)) dt - NA^{rj) 


'0 i=l 


4a 




Y1 


dxf{s). 


Note that we have |Pe(a;)| < i and therefore we can deduce 


lim sup lim sup lim sup 

7j—>-0 e—>-0 N^oo 

^ log Eat exp dt - NA^{t]'^ } ~ ^ 

for all a > 0 by the exactly identical way to the previous step. Thus, we can conclude (I2.46P as 
well. □ 


3. Exponential Tightness 

In this section, we establish the exponential tightness of which can be deduced from 

the following result. 


Theorem 3.1. For any e, a > 0, 


lim sup lim sup — log Fn 

( 5^0 N^OO 


✓ ' 



i : sup \xf{t) — > e 

> 

> Na 

0<s,t<T 



. h-i|<5 




— 00 . 


Before proving Theorem 13.11 we briefly explain the reason for the exponential tightness of 
being a corollary of this theorem. We can prove exponential tightness of 

by showing 


lim lim sup — log Fn 

5^0 Af—>00 


sup 

0<s,t<T 
|s—1|<(5 


1 

N 


(J=(xf (t)) - W)) 

C= 1 i&lN 


> e 


= —00 
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for any (Ji, J 2 , • • • , Jm) £ C(T)™' and e > 0. It is obvious that this estimate is a direct consequence 
of 


1 


lim lim sup — log Pat 

(5^0 N-i-oo -/V 

|s—1|<<5 

for each c. We can deduce this estimate from Theorem O since 


0<s,t<TtV 


i£l^ 


= —00 


Pa 




0<s,t<T, |s-t|<<5 

< Pa 


i: sup |a:f(i)-a:f(s)| > - 

0<s, t<T, |s—i|<(5 ^ 


> 


Ne 


We return now to Theorem 13.11 The basic strategy is to divide the estimate into the normal and 
small time regimes as before. To carry this out, we first observe that 


where 


i: sup |a;f(t)-Tf(s)| > e I C 7?]) US’|^ 5 ([r/, T]) 

0<s,t<T, |s-t|<(5 J 


SeAlv, T]) = <i: sup \xf{t) - xf (s)| > e 

y ? 7 <s, t<T, |s—1|<(5 

'S'e, 5 ([ 0 , r]]) = I'i- sup \xf{t) - xf (s)| > e 

I 0<S, t<77, |5 —1|<(5 


(3.1) 


(3.2) 


for all T/ > 0. Consequently, Theorem 13.11 can be separated into the following propositions. The first 
one is the normal time regime type of estimate. 


Proposition 3.2 (Normal time regime). For any rj, e, a > 0 

lim sup lim sup — log Pat [|S'e^ 5 ([ 7 /, T])| > Na] = — 00 . 
5—>-0 TV—>-oo 

For the small time regime, we have 

Se,si[^^ r/]) C If : sup |xf(t)-xf(0)| > ^| 

[ o<t<v ^ ] 

and hence it is enough to prove the following estimate. 

Proposition 3.3 (Small time regime). For any e, a > 0 


(3.3) 


lim sup lim sup — log Pat 
?7^0 TV—>-oo ^ 


i : sup |xf(t) -xf( 0 )| > 

0<t<ri 


> Na 


= — 00 . 


For the normal time regime, it is possible to transfer the estimate to that of the equilibrium 
process bv usiner Lemma 12.141 Then, we can applv the well known methodoloerv {e.g., m) based 
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on Garsia-Rumsey-Rodemich’s inequality to the equilibrium estimate by making a small adjustment. 
However, for the small time regime, we cannot send the estimate to the equilibrium and therefore 
have to adopt a different approach. 


3.1. Normal Time Regime. By setting f{x) = xi in (II.2p and (II.3p . we obtain 


xf{t) = Pi{t) + Af{t). 


(3.4) 


where 

Af{()= (3.5) 

which can be regarded as the difference between the left and right collision times for particle xf{-) 
and also measures the deviation of the lifted particle xf{t) from the underlying Brownian motion 
/3i(t). In contrast to the averaged local times in (12.ip . (12.2p or (j2.3p . the behavior of Af [t) is 
unacceptably noisy. Thus, we now present a way to control this object. 

By dHT}, the estimate (13.3p can be divided into 


lim sup lim sup 
5^0 N^oo 


^logPTV 



sup \l3i{t) - f3i{s)\ 

r)<s,t<T, |s—1|<5 



lim sup lim sup 
(5^0 N^oo 


> Na 


= —oo. 


^logPiv 


< i : sup 


-1 

> Na 

? 7 <s, t<T, |s—1|<5 


i 



—oo. 


(3.6) 


(3.7) 


First of all, (13.6p is standard because /3i’s are independent Brownian motion. The main challenge 
is (j3.7p . Since the event inside the bracket of (|3.7p is invariant under the permutation of labels, we 
can apply Lemma [2.141 to send the estimate to the equilibrium as following proposition. 


Proposition 3.4. For any e, a > 0 and T > 1, 


lim sup lim sup — log 
(5^0 N^oo ^ 



sup 

0<s,t<T 
|s—1|<5 





> e 


> Na 


—oo. 


Remark 3.5. We expanded the time window from [r/, T] to [0, T] to reduce unnecessary notational 
complexity. We also assumed T > 1 without loss of generality. 


3.1.1. Garcia-Rumsey-Rodemich’s Inequality. For (j) & C'([0, T], M), let us define 


5t((/.) 


sup sup 

0<5<i 0<s,t<T 
|s—t|<(S 


\(j){t) - (/)(s)| 
^logi 


then we have the following result. 
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T rT 


Lemma 3.6. For T >1, we have 

5r(0) <Ci + C2log r [ 
Jo Jo 

for some positive constants Ci, C 2 - 


exp 


(fit) - (fis) 




dtds 


Proof. Define p(x) = , ^'(x) = — 1 and M = exp | | dtds so that M > > 1. 

For |t — s| < S, by Garcia-Rumsey-Rodemich’s inequality (c/. Section 1.3 of [36]), 

r\t-s\ r 4 (M - T^) 1 
\(j)it) - (fis)\ <8 log<l-F- ^ - >dp{u) 

-3 f 4M\ , 

< 2 y M 4 log ( M -I- ^ \ du 

= 854 logM + du. 

Therefore, the proof is completed since we have 

f 


j u 4 log -I-drt < [ u 4 (^2 + 21og —^ du = 4054-|-854 log 


□ 


3.1.2. Proof of Provosition [2fTR We return now to Proposition 12.12] For 5 < i, we have 




<P^« 


< i: sup 

A«{t) - Af (s) 

-1 

> Na 

i 0<s,t<r, |s-t|<(5 


i 



N 


i=l 0<s,t<T, |s-t|<5 


<^7V 


N 


i=l 


EsRA" > 


Nae 


> Nae 


^logT 


< exp < — 


Nae 


Cs^logl 




r , N 


i=l 


where C 2 is the constant from Lemma 13.61 Therefore it suffices to show 

N 


lim sup ^ log 
Af^oo tv 




2=1 


< C 


(3.8) 
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where C is a constant which does not depend on S. By Lemma l3.6[ 


(A^)| 

nil 


<IE^' 


Liiv 






dtds 


= eC'2 


£ijv 


Cl rT rT 

hiAT / 


P I " Uf ((i) - Af {.i) 


dtidsi ■ ■ ■ dt^dsN 


S~^ />T ^eq f (ti) — A^ (si) . 

V / ••• / E^exp^Vei--^ }> dti • •-dsAT. 

c,=±i,w^o io ' 


. 2=1 


^ti - Si 


We will prove the following lemma in the next subsection. 


Lemma 3.7. For any ai and 0 < Si < ti < T, 




exp I'^at (^Af^{ti) - Af (si)) 


. 2 = 1 


N 


< exp + af){ti - Si) 


2=1 


where C is a constant only depending on T. 

By assuming this lemma, we can bound (13.91) by 


5 C 2 


rT fT rT ( ^ 

E / ■ p-p cE 

.±ivPo Jo Jo [ 


{Vti — Si + 1 ) 7 dtidsi ■ ■ ■ dtNdSN 


and hence (|3.8p is proven. 

3.1.3. Proof of Lemma \d. 71 The final step to prove Proposition 13.41 is Lemma 13.71 We prove this 
lemma by a series of estimates. 

Lemma 3.8. For any ai and 0 < s* < L < T, 


(3.9) 


(3.10) 


exp < oji (^Af{ti) - Af(si)j - ^ {Af(ti) - Af{si)) 
I 2 = 1 *- 

Proof. Let us define V{t, x) = «i](^)^*(^) where 


< 1 . 


(3.11) 


V^{x) = Y 

hi 7 ^* 


aj 


ai (<5+(xi - Xj) - 6+{xj - Xi)) - {6+{xi - xj) + 6+{xj - Xi)) 




















LARGE DEVIATION PRINCIPLE FOR INTERACTING BROWNIAN MOTIONS 


28 


Note that we can rewrite ()3.1ip as 

E^exp|y V{t, < 1. 

Now, as in the proof of Proposition 12.151 we can obtain 
lE^^exp |y V (t, x^(t)) dtj 


(3.12) 


(3.13) 


< exp 


ir 


sup 


V(t, x)f(x)dx - Piv(/) 


Si fe^jvnv(^jv) UGjv 

by Feynman-Kac’s formula and the variational formula for the largest eigenvalue of + V. Note 
that we can bound V(t, x)f{x)dx — T>fq{f) by 

a? XN / r— r— \2l 

\fij ~ fji\ ~ ^ 


E 


i¥=j ^ 




{x)dSij{x). 


It is not difficult to check the last expression is non-positive because of the elementary inequality 

1 „ / _ 2 


x'^{a + b) +-y'^ (^y/a — VbJ >xy\a — b\ 
for a, 6 > 0. Thus, the RHS of (13.131) is bounded by 1 and hence ()3.12p holds. 


□ 


The next estimate is a stronger version of Lemma 12.51 


Lemma 3.9. For any 1 < i < N and f G we have 



Proof. We introduce a function a^^\x) on Gm for k ^ i hy 

N 

^k\x) = \0,xu-Xi]{Xj - Xi) 

i=i 

which counts the number of particles between x* and Xk in the clockwise sense. We remark here 
that this function also appeared in |10| to estimate Af{t). We normalize a^^\x) by c^*^(x) = 
(7^^\x) — and set cf\x) = 0 for simplicity. We can define a piecewise constant vector held 
Cj(x) = and apply Green’s formula (12.5p . First note that V • Cj = 0 and thus we 

only need to concern about boundary terms. On F^i with k, I ^ i, we have 

(Ci(x), Ck - ei) = 4*^(x) - c4(x) = -1 

because x; = -|- 0 on Fi^i and hence 4*^(x) = (T^^\x) + 1. For the boundary Fji, 


N-2 
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because = iV on Fji and {x) = 0. Similarly, on Fij, 


(Ci(x), Ci - Cj) = cf'\x) - c'-\x) = 




N-2 


- 3 ! V-/ 2 

We now apply Green’s formula with the vector field Cj(x): 

N 


^ [ [Vkf{x)]ckix)dx 

k=iJGN 

= - ^ f fki{x)dSki{x) 


IV- 2 

3-3T^ 


fij{x)dSij{x) + / fji{x)dSji{x) 


f ifuv + fvu){x)dSuv{x) + ^ [ ifij + fji){x)dSij{x). 


Since we have bound |cfc(x)| < —we can derive 

X] [ ifij + f3i)i^)dSijix) 




N 

ifuv +fvu)ix)dSuvix)+ '^ f \Vkfix)\dx. 


N - JF 


_i JGn 


k=l 


Note here that the RHS is bounded by 2N + 2y^8A^PAr(/) due to Lemmas 12.41 and 12.51 


□ 


Lemma 3.10. For any a* > 0 and 0 < Si < U < T 


N 


N 




. i=l 


. i=l 


Proof. It suffices to prove 


'Si 


exp < / Vi{x^{t))dt > < exp {8A^(a + a^)(L — Sj)} 


(3.14) 


for any i and a > 0, where Vi{x) = aYlj-j^^i l^+i^i ~ Xj) + ^+ixj ~ ^J*)] • Note that the LHS of 
(j3.14l) is bounded above by 


exp I {ti - Si) sup S “ y] / (/l(®) + fji{x))dSij{x) - VnU) 

' fe^Nnn.^N) I j.jjLiJFij 
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as before. Finally, by Lemma 13.91 

“ X] / + fji{x))dSij{x) - VN{f) 

<a{2N + V32iVP^) - V^if) 

< SN{a + a^) 

and we are done. □ 


Consequently, we are able to prove Lemma 13.71 by Lemmas 13.81 13.101 and Cauchy-Schwarz’s 
inequality. 


3.2. Small Time Regime. In this subsection, we provide a detailed proof of Proposition 13.31 The 
small time regime differs from the normal time regime in that Af {t) on t G [0, rj] cannot be properly 
controlled if the process does not start from the neighborhood of the equilibrium. Therefore, it is 
not possible to work directly with xf{-) as in the normal time regime. Instead, we introduce an 
intermediate process zf^{■) where 

Zi{t) = xf{t) + ^ ; l<i<N (3.15) 

where p(x) = x on [0, 1]. These adjusted processes were introduced in |10j and were turned out to 
be martingales with respect to the same filtration with x^{t). More precisely, we can prove that 

zf{t) - zf{0) = ^f{t) + ; 1 < z < iV (3.16) 

where 

k:k^i 

For the details, see Proposition 2 of m- 

We first develop the exponential tightness of {z^ as an intermediate step, which can be 
formulated as follows. 

Proposition 3.11. For any e, a > 0, 

lim sup lim sup — log P tv 

rj^O N—>-oo ^ 


sup 

0<t<r) 


|2f (*) - if (0)1 > 


> Na 


= —oo. 
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Proof. By (I3.16p . it is enough to show that 
lim sup lim sup — log P n 

ri^O N—>-oo ^ 


lim sup lim sup — log Ptv 
r]^0 N—>-oo 


< i : sup 

Pfit) 

y 0<t<r] 



< i: sup 


y 0<t<ri 



> e 


> Ne 


> Na 


> Na 


= —oo 


= —oo 


(3.17) 

(3.18) 


respectively. Let us first consider the (13.171) . We can easily bound this as 


Pat 

< i : sup 


-1 

> Na 

< Pa 

■ N 

sup \/3i{t)\ > Nae 


y 0<t<ri 


J 



0<t<r, 


Since {/3i{t)}^i are independent Brownian motions, it is easy to check that this probability is 
super-exponentially small. 

The next step is (|3.18p . We define two adapted processes 

and then we can rewrite (|3.18p as 


lim sup lim sup — log P n 

ri^O N^oo N 


N 


Y,Ctiv)>Na 


. 1=1 


= —oc. 


By Chebyshev’s inequality it suffices to show that 


1 f ^ 1 

lim sup lim sup — log Eat exp < a / Cf {rj) \ <C 
77^0 N^oo N y ^ J 


(3.19) 


(3.20) 


for all a > 0 where C is a constant does not depend on a. We shall prove (I3.2nh with C = log 2. 
First we prove 

1 n 

(3.21) 


(-+{„) < J- 


f\-{s)dM^{s) . 

Jo 

To see this, we only need to concern about the case where Ci'(^) = 1- Define 

rjo = inf |t : {f) > A^ej 


and then r]Q < r/. Thus, 


__ rvo _ 

/ C.-{s)dM^{s) = / l-dMf{s) 

Jo Jo 


Mtido) 


> Ne 
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due to the right-continuity of the jump process and we proved (|3.2ip . We return now to (|3.20l) . By 

(HSU), 


E7vexp|a^C»^(r/)| <IEivexp|^^ ^ C- 


-is)dMfis) 


< 


ei=±l,Vi 

Note that we can rearrange each summand in (j3.22p in a way that 


Y (s)| . (3.22) 

=±l.Vi i=l "'O J 


Eat exp <1^^ I tiCi is)dMf^{s) 


N 


2=1 


ETvexpJ^ Y [ ^iCi {s)d[Mjl{s) - M^i{s)] 


En exp i Y [ 
1 l<i^k<N''^ 




Ne 


-dMiUs) 


(3.23) 


where Uik{s) = eiCf{s) — ekCk(s). Then, since each M^{t) is the compensated Poison process with 
rate \NA^f^{t), 

N 


Etv exp Y^ 
i,k=l 


/ 

Jo 




™ 2a^ _ 20 ^ _ \NdAS.{,) 


= 1 . 


Therefore, by Cauchy-Schwarz’s inequality, (|3.23p is bounded by 


Eat 


N 


expY exp^^^-- — -1 ) AA^d24i;i,(s) 


2,fc = l ' 


Ne 


Ne 


1 
2 
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For sufficiently large N, we can bound the last expression by 


En 


N 


< Eat 


= Eisf 


exp / I 

N 

exp ^ f 


f 2auik{s) 


Ne 


XNdAfiis) 


^ 16a^A , ,Ar 


i,k=l' 

Iba^A 

exp—^ 


iVe2 


dAg{s) 


NA^{r]) 


since |ujfc(s)| < 2. These series of estimates enable us to bound 

f 1 r c 1 r^2 

Eat exp a ^ (f iv) [• < 2^EAr 


2 = 1 


IGn A AT, , 
exp <; — 2 — NA^{tj) 


Now the proof is completed by Proposition 12.171 


□ 


Now we prove the tightness of by starting from that of The methodology 

for this step was developed in Proposition 3 of m for a fixed i. The situation here is slightly different 
but we can still burrow the core idea. 


Proof of Proposition \3.A We start by defining two stopping times and r,- ^ as 


for each i and then 


I 0 <t<r; 

Thus, it suffices to show 


= inf {t : xf{t) - xf (0) > e} 

= inf {t ■ xf{t) - xf (0) < -e} 


sup |xf (t) - xf (0)| > 4 = < ?/} U |r._^ < r/| 


1 


lim sup lim sup — log Pat 

27—>-0 N—^OO ^ 


{i ■■ r+ < r?} 


> Na 


= —00 


since r- can be handled by the exactly same manner. Let us dehne 


(3.24) 


N 


U^it) = 


1V(A + 1)^ 


^p(xf(t) -xf(t)) 
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SO that {t) = xf{t) + uf {t). Then, 

= {<e < 


c 


<«e}u{r+<??, zf(r+)-zf(0) > Ke} 
{r+ < r/, |e + ixf (r+J - uf (0)| < ne] U | ^sup^ \zfit) - zf (0)| > Aej 


for any k > 0. Note that the second set is super-exponentially negligible by Proposition 13.111 For 
the first set, we have 

e + uf {t+^) - {0) <Ke 




C 


< r], 


} 


{r+ < V, uf{T+J - nf (0) <{k- l)e} . 
Take S (^““(T) satisfying ^{x) < (/>e(x) < p(x — e) + (1 + K)e so thal@ 
nf(r+)-irf(0) 

= jv(A + l) ^ ^ (<J) - ^ (®f (0) - (0)) 

^ i=i 


JV(A+l)^ 


N 


> 


iV(A + l)^ 


^ U, ix^irt^,) - xf(0)) - 4>e (xf (0) - xf (0)) - (1 + K)t 


Thus, — uf^ (0) < (k — l)e implies 


N 


J^(^X + 1) ^ “ xf{0)) - (/>e(xf (0) - xf (0)) 


/ ^ K + 1 7 

< (k — l)e + --e := —--e. 

“ ^ ^ A + 1 A + 1 

We choose k small enough so that 7 > 0. Then, the RHS of (13.251) is a subset of 

N 


sup 

0<t<ri 


[(Pe {xfit) - xf (0)) - (j)^ (xf (0) - xf (0))] 


i=i 


> 76 


(3.25) 


(3.26) 


^The existence of such a function is proved in Proposition 4 of [10] 
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By Ito’s formula, 


N 


N 


[</>6 - xf (0)) - (j)e {xf (0) - xf (0))] 


i=i 


< 


N rt 


- xf (0)) ds 


i=i 


+ 


N rt 




<!l 


+ 


N rt 


^Y [ (®f (0)) dPj{s) 




for t < rj. Note that the local time does not appear since the expression (I3.26P is totally symmetric 
with the function — xf (0)). Thus, (13.261) is a subset of 


sup 

0<t<r) 


1 ^ 

^Y^J^ “ a:f (0)) d/3j(s) 


>11Vn 


(3.27) 


for sufficiently small r/. We now regard Yj=i fo (®) “ ^f (0)) d/3j{s) as a time change of 

Brownian motion 

- xf^{0)fds 

where Bi(-) is a Brownian motion starting from 0 under Pjy Since 




s) - Xi (0)) ds < 


T] := CeV 


the event (j3.27p is a subset of |supo<t<c'^^ \ 
Finally, it suffices to show 


1 


lim sup lim sup — log Pat 
ri^O N—^oo ^ 


> Na 

to complete the proof. However, this is obvious since we have a trivial bound 


i : sup \Bi{t)\>—VN 

0<t<C^ri 2 


= —oo 


Pat 


i: sup \Bi{t)\>^VN 
0 <t<CeV ^ 


> Na 


N 

2=1 


sup |Hj(t)| > ^Vn 
0<t<c^r] 2 
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for large enough N, and then by a property of the Brownian motion 
Pat 


sup 

0<t<CeV ^ 


8^/C^ 

< - , exp 


7e\/27riV 


2 V JV I 


2C,rj 


□ 


4. Diffusion of Colors 


4.1. Introduction. In this section, we develop the LDP for the empirical density of colors. We 
recall from Section [1.2.21 that the empirical density for colors {jl^ (t) : 0 < t < T} is defined as 
(11.101) which can be regarded as a Markov process Qjv on T], First, we state the 

hydrodynamical limit theory of 

Theorem 4.1. Suppose that Assumption 2 is satisfied with the uncolored initial measure p^{x)dx for 
some bounded Then converges weakly to Qoo which is a Dirac mass concentrating 

on the single-valued trajectory 

{p(t, x)dx = {pi{t, x)dx, P 2 {t, x)dx, • • • , pm{t, x)dx)^ : 0 < t < T} 
where p{t, x) is the unique weak solution of the partial differential equation 

^='^V-[D{p)Vp] (4.1) 

with initial condition fp{x)dx where the diffusion matrix D{p) is given by 

Remark 4.2. An extension of this hydrodynamic limit result is obtained in |30| which consider a 
system of two mechanically different types of particles. 


As we already observed in (|1.13p and (|1.14l) . each component of equation (|4.ip can be written as 


dtpc = 2 ^ 


^ ^ Vp ■ 

T—^Pc + T—Pc 

A + p A + p 


; c=l, 2, 


m. 


(4.2) 


where p = Pc- Since p is the solution of the heat equation dtp = ^Ap with the initial condition 
p^{dx) = ^4.21) is just a usual linear parabolic equation with smooth coefficients. 

Thus, the uniqueness of the solution is immediat^ and consequently Theorem 14.11 is a corollary of 
Theorem 14.181 The final comment regarding Theorem 4.1 is that we did not need assumptions on 
f^{dx) more than what is stated in Assumption 2 because our estimates in Sections 2 and 3 did 
not impose any further conditions due to our careful analysis on the small time regime. 

The next step is the large deviation theory for {QAf}^=i- In order to concentrate on the large 
deviation of dynamic evolution and simplify the arguments regarding the initial deviation, we as¬ 
sume the following throughout Sections 4 and 5 in addition to Assumption 2. 


®This uniqueness is also a direct consequence of Theorem 14.201 
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Assumption 3. The initial configuration of particles is i.i.d. with a bounded probability density 
function p^{x) on T. 


Remark 4.3. Reading our proof carefully reveals that the LDP is still valid under many general 
initial configurations, e.g., deterministic configuration. However, the boundedness assumption on 
p^{x) is essential, especially when we establish the compactness property of the rate function in 
Lemma 14.91 


We now state the LDP for {QAr}^=i under Assumptions 2 and 3. First of all, for each color c, 
C ^(T) is defined by 


^°(T) = |r(x)(ix : J r{x)dx = pc, r{x) > 0 

where p^ is given by (|1.8I) . Then we will show in Section 4.2 that the LDP rate is infinite outside 
^([0, T], where 

m 

■^co/or-(T) = n-^°(T). 

C=1 

In the domain (^([O, T], -y^color^))■: function is = Ifffifi-) + where the 

dynamic rate function ^Tyni-) is defined by 


^dyn 


(p(-> 



^V-[D{p)Vp] 


2 

dt. 

-hA{p) 


(4.3) 


for p{-, x)dx G (^([0, T], -y^color^)) where A{p) is defined by (ll.lGjl . Here, it is necessary to clarify 
exactly what is meant by the RHS of (|4.3I) . The H-i a norm can be explained by the variational 
formula (c/. (2.24) of [2i] f 


sup < / 4>^p{T, x)dx — / <))i/5(0, x)da 

</> Jt Jt 


(4.4) 


+ 


'0 Jt 


p)(k\ 1 1 


{t, x)dxdt 


where the supremum is taken over (j) G C°°([0, T] x T™). However, the last expression is still not 
well-defined as it stands because it involves Vp, which might not exist. Thus, our starting point 
should be to obtain a reasonable explanation of (|4.4I) and the basic properties of this rate function, 
e.g., the compactness and lower semicontinuity. Section [4.21 is devoted to this project. Then, we 
will establish the large deviation upper and lower bounds in Sections 14.31 and 14.41 respectively. 


4.2. Rate Function. 

4.2.1. Well-definedness of Variational Formula To define the rate function in the sense of 

(14.4p . we need some a priori regularity result for p(-, •), as well as energy estimates on the domain 
of the rate function. Recall here that p = (pi, p 2 , ■ ■ ■ , Pm)^ is the m-dimensional vector of the 
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density of colors and p = Pc denotes the total density. To begin with, let us define a set 

^color C C([0, T] , -^color^)) consists of p{t, x)dx which is weakly differentiable in x with the 
energy estimate 

/ (yX^xip)'^p) (t, x)dxdt < oo 
and satisfies the finite initial entropy condition 


/ p(0, x)logp(0, x)d: 
Jf 


.X < oo. 


(4.5) 

(4.6) 


Note here that we can compute Vp^x^x(/o)Vp explicitly as 


^ + P ^ (-^ + P)Pc 

and accordingly, (14.51) is equivalent to 

i{p) = Jq fj (t, x)dxdt < oo and, 

/c(p) = /t (S)1 ^)dxdt < oo 1 < c < 


(4.7) 


m. 


In particular, the finiteness of I(p) implied by (14.51) implies that p G T2([0, T] x T) or equivalently 
p € T2([0, T] X T™) as follows. 

Lemma 4.4. Suppose that p is weakly differentiable and satisfies I(p) < oo. Then p G ^ 2 ( 6 , T, Loo(T)) 
and the L 2 norm is bounded by 2{I{p) +T). In particular, p{-, x)dx G implies that p G 

L2([0, T] X T™). 

Proof. Let be the heat kernel on T at time and pg = p* fie- Then, 


Pe{t, x) - Pe(t, y) = [ Vpe(t, z)dz < \ [ (t, z)dz 

J[x,y] y Jt 

by Cauchy-Schwarz’s inequality. Integrating the above expression against y give us 


sup p, 


ft, x) < ^-(t, z)dz + 1 < ^ 


{Vpf 


{t, z)dz + 1 


This implies sup,j.gj p^(t, x) < 2(I(p) + T) and we can obtain the desired result by taking 


0 . 


□ 


Remark 4.5. Henceforth, for any function / on T, fe denotes f * fie where fie is the heat kernel on 
T at time e^. 


The next two lemmas prove that the domain of the rate function is included in 
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Lemma 4.6. For any G E T] x T), 

rT ^ 


limsup — loglEAT exp < f xf (t)) — xf (t))^ dt\ < 0 (4-8) 

N—^oo -A I Jq I 


and for each color c, 

lim sup lim sup 

V—^oo 


^logEArexpj^ ^ VG(L xf{t)) - 12G^(t, xf{t)) ^1 + 


< m 


ieiF 


dt 


(4.9) 


Proof. For (j4.8p , the expression is symmetric and thus we can consider the model as non-interacting 
case. Then, 


log Eat 
N 


(t)) - 2G^{t, xf (t))] dt 

^logEAT exp I f [VG{t, xf{t)) - 2G‘^{t, xf{t))] dt 
i=i L Uo 


(4.10) 


If /3 is a standard Brownian motion under P, then by Feynman-Kac formula and the variational 
formula for the largest eigenvalue, 

r-T 


log E 


exp 


[ [VG{t,(3{t))-2G\t,P{t))]dt 
Jo 


< 


[ sup < [ h{y){VG-2G^){t, y)dt 

Jo h£C°°{T),Lh{y)dy=l,h>0 I ./T « ^\V) I 




0 hec°°(T),/.,, 

T 


sup 

0 hGC°° (T), fj h(y)dy=l, h>0 




4G(t, y) + 


^hjy) 

Kv) 


dy > dt 


< 0 . 


Therefore, (14.101) is non-positive. 

For (14.91) . we should take the interaction into account and this requires us to consider additional 
Pi,e{x^{t)) part. In the spirit of Lemma 12.141 we can replace Etv in (14.91) by Then, the 

estimate with E^^°^ is equivalent to the one with E^ since 



< 

cUPty 


dr^ 

^color 

CXD 

pi^color 

CO 

dP^v 


< rn^ IIpoII 


N 

oo 
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by Assumption 3. Consequently, we can substitute E^r in (14.91) by E^ and the price of this substi¬ 
tution is m 11 Pol loo- equilibrium estimate, by the standard argument as before. 


^logE^'expl / 


VG{t, xf(t))-12G^(t, xf(t)) 


2/-+ ^N/ 


s r / /( 

Jo fe^N [ Jgn 


^ ^ Pi,eix 

PiA^)) 


dt 


VG{t, Xi) - UGAt, Ti) 1 + 


A 


dx 


-Mill 

N I 


Now we apply Green’s formula (|2.5p with the vector field Yliei^ G{t, Xi)ei such that 

V [ f{x)VG{t,Xi)dx = Ui + U 2 

TTmJGm 


iei^ 


where 


Ui = 


f Vif{x)G{t, Xi)dx 

\i-TN •’Gn 




U2 = EE/ G{t, Xi){fki{x) - fik{x))dSik{x). 

i£lN k-.k^A Gk 


We can estimate Ui as 


c.<E/ [ 


iei^ 
< 12 


12f{x)G‘^{t, Xi) + 


1 (Vi/(x))^ 


48 fix) 

V [ f{x)G‘^it,Xi)dx + ^VNif). 

Arir Jg^ 6 


dx 


ie/A 


For U 2 , we can apply Proposition 12.11] such that 


(4.11) 


(4.12) 


r o _ _ \i\j 

+ ^/J^fiAG'^it, Xi) + —{y/JA-y/fkAix)dSikix) 

^ [ ifik +fki)ix)G‘^it, Xi)dSikix) + ^'DNif) 

iei^,k^i'^Gk 

< y ^ fix)peAAG‘^it, Xi)dx + CNe^ ^1 -h ^ + ^'^nU) 


(4.13) 
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where C is a constant only depends on G. Thus, we have by (|4.12p and (14.131) that 



Thus, the proof is completed. 


□ 


Based on the previous lemma, we can restrict the domain of rate function to Note that 

{Qiv}?/=i satisfies the LDP due to Bryc’s inverse Varadhan Lemma(c/. Theorem 4.4.2 of [3]) because 
of the exponential tightness result of Section 3. Hence, we will temporarily denote the rate function 
by /"*(•) in the next lemma, since we do not know exact form of the rate function at this stage. 


Lemma 4.7. If < oo, then fi{-) = {fji, /U 2 , • • • , G ^JSior- 


Proof. First of all, we know from the Theorem 3.1 of |18) that < oo only if fi{t) is abso¬ 

lute continuous with respect to the Lebesgue measure for all t. Let us define a functional He on 
C([0, T], .^(T)) by 

EcifLi-)) = [ dt [ [VG(t, x) - x)] //(t, dx) (4.14) 

Jo Jt 

where G G T] x T) and ^ = Y1T=1 d-c- Then (14.8p can be rewritten as 

lim sup ^ log[exp {A^HG(/i(-))}] < 0 

N^oo N 


and therefore by Varadhan’s Lemma, He(/1(-)) < particular, if /!(•) = p{-, x)dx satisfies 

< oo, then we have 

f f VG(t, x)p{t, x)dxdt < + 2 / f G‘^{t, x)p{t, x)dxdt. 

Jo Jt Jo Jt 


and therefore 


f f VG(t, x)p{t, x)dxdt < gJ f f G‘^{t, x)p{t, x)dxdt (4-15) 

Jo Jt V 4o Jt 

for some constant G > 0 which does not depend on G. If we define the inner product (•, •) on 


G([0, T] X T) by 


(Fi, F 2 ) = [ [ FiF2pdxdt 

Jo Jt 

and let Lp([0, T] xT) be the Hilbert space by taking completion and equivalent class. Then (14.151) im¬ 
plies that the functional /(G) = fj, pVGdxdt is a bounded linear functional on G^4([[o, T] x T) C 
Lp([0, T] X T). By Hahn-Banach’s theorem, we can extend /(•) to and then Riesz representation 
theorem gives us a function H £ such that /(G) = (G, H)^ . Therefore, p is weakly differentiable 
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with V/? := —Hp. Moreover, since H E T] x T) we obtain 

f H^{t, x)p{t, x)dxdt = f f dxdt < oo. 

T Jo Jt P 

This proves the finiteness of I{p). 

Likewise, we can derive from (|4.9I) that 

limsupy J (VG - 6G‘^ (^1 + pcdxdt < I'^ip) + m\\po\\oo- 

Since p E L2([0, T] x T) by the finiteness of i{p) and Lemma (14.41) . the LHS is 

and the RHS is independent with G. Therefore, we can repeat the previous argument to prove that 
Pc is weakly differentiable and Ic{p) < oo. 

Finally, if < oo then by Sanov’s theorem fj. p(0, x) log ^-^f^dx < oo and this implies the 

finiteness of the entropy (I4.6jl . This finishes the proof. □ 



Now, we are ready to explain the dynamic rate function in the sense of (|4.4p for p{t, x)dx E ^^lor- 
The only part which is not well-defined in (14.4p is 

rT 


[ f V(tJD{p)Vpdxdt. 

Jo Jf 


First observe that the cth element of D{p)Vp is y^Vpc + and hence it is enough to show 

the finiteness of 


II 


^ Pc 


X + p 


dxdt and 


II 


Pc 


X + p 


Vp 


dxdt 


for p(t, x)dx E The first one is bounded by 


f I 


(Vpc 


T (A + p)p, 


-dxdt 


f I 


Pc 


X + p 


dxdt 


and therefore finite for p(-, x)dx E ^^icyr- The second one is bounded 


1. 

2 


by 




PcP 

(A + p)2 


dxdt 


1 

2 


which is also finite since — P- Therefore, for p(-) =p(-, x)dx E can define 

^dyniPi'i x)dx) through the variational formula (14.4p . We finally set dTyniK-)) =ooforp(-) ^ ^ color' 


Remark 4.8. ffenceforth, we write p E and /^„(p), instead of p(-, x)dx E 

/^^(p(-, x)dx), respectively, for simplicity. 


4.2.2. Lower Semicontinuity. The next step is to establish the lower semicontinuity of the functional 
dcolori') equivalently To carry this out, we start from a compactness result. 
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Lemma 4.9. Suppose that p G O'^d < oo. Then, we have 

f f '^P^xAxip)'^pdxdt<C{l + I^i^^{p)) (4.16) 

Jo Jt 

for some constant C. 


Proof. For g G (^““([O, T] x T™) with /.j. x)dx = 0 for all t G [0, T], we can consider a semi¬ 
norm 11511^ i(A(p)) ~ fo Ilfl'II-i A{p) taking completion and equivalence class, we obtain 

space. For h G (^““([O, T] x T”^), we have another semi-norm = JjVh^ A{p)'Vhdxdt 

and we can obtain space in a similar manner. These two spaces are dual each other and 
hence, for g G J^-\ and h G the integral x)h{t, x)dxdt is well-defined and satisfies 

Cauchy-Schwarz’s inequality 

/ / 9{i, x)h{t, x)dxdt < \\g\\.^_^(A(p))\\h\\.^i{A{p))- 

Jo Jt 

Now, we write 


log/5 = (logPi, logp 2 , • • • , logpm)^ 

G = dtp-]^\/[Ax{p)yp\ 

then G, V [^x(p)Vp] G J^-\ and logp G where the norms can be easily computed by the 
variational formula given in (2.24) of |24j such that 


\\G\\l,_,^APp))=l7yn{p) 

(4.17) 

V [Ax{p)Vp\ \\Y{A{p)) = J^Vp^xAxip)'^pdxdt 

(4.18) 

IUogp||^i(A(p)) = J^Vp^xAxip)'^Mxdt. 

(4.19) 


Now, let us consider the entropy functional 

m 

Ht{p) = Y. 

c=l 


Pcit, x)logpcit, x)dx 
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for p E Then, by (|4.17p and (14.191) . 

Hrip) - Hq{p) = (logp)^ Qv • [Axip)Vp] + dxdt 

^ P + f f {log p)"^ Gdxdt 
^ Jo Jt Jo Jt 

f [ +IUog^ll^i(A(p))l|G'||^_i(A(p)) 

^ Jo Jt 

-~\fo 

Note that 

Ho{p) < f p{0, x)logp{0, x)dx < [ p(0, x)log^^^^y^(ix + log||/)o||oo 
Jt Jt Po{x) 

and Jrj, p(0, x) log dx is the large deviation rate for /U'^(O) and therefore bounded by I^^ip) 
by the contraction principle. This proves (14.16^ . □ 


Establishing the lower semicontinuity or I^ig^{-) requires a few convergence results. The following 
lemma implies that the weak convergence can be combined with some energy estimates to obtain 
the strong convergence. This lemma is motivated by Lemma 4.2 of [24| but our formulation and 
proof are differ slightly in that we do not have a priori boundedness of the density. 


Lemma 4.10. Suppose that {fN{',x)dx}^^^ C (^([0, T], ^{T)) satisfies 

fN{', x)dx /(•, x)dx weakly in ^([O, T], ^{T)) 

/ /Ar(L x)dx = / f{t, x)dx = / for all t, N 
Jt Jt 


^ '-{VfN? 


f I 


n onIn 

where positive functions ai\f{-, ■) satisfies 


dxdt < G for all N 


/ aN{t, x) 

Jt 


< M 


(4.20) 

(4.21) 

(4.22) 

(4.23) 


uniformly in t, N for some M > 0. Then, fn ^ f strongly in Li([0, T] x T). Moreover, if 
fNi f £ T2([0, T] X T) and = 1 for all N then /jv —>■ / strongly in L2([0, T] x T). 


Proof. For the first part, we first recall the notation of Remark 14.51 and then it suffices to show 

limlimsup [ [ \ fN - (/a)J + \ifN)^ - fe\ + \ fe - f\ dxdt = 0. 

e^o AT^oo Jo Jt 
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For the second term in this limit, observe first that (/Ar)^ fe pointwise as —>■ oo due to the 

weak convergence and 


[ [\{fN)^\dxdt= [ [ \ fe\dxdt = Tf. 

Jo Jt Jo Jt 


Therefore we can apply Scheffe’s Theorem to check the desired convergence. The third term obvi¬ 
ously tends to 0 as e —)• 0 and therefore it is enough to show 

rT 


lim sup / 
^^0 N Jo 


[ l/v - (/v), 

Jt 


dxdt = 0. 


By Cauchy-Schwarz’s inequality we can bound 

r-T 


, , < 
Je\ — 


[ [ \ fN - {In) 

Jo Jt 

Now we can bound fj {^/Jn — ^ dx by 

fwit, x) - \/fNit, x + y)^ (j)eiy)dydx 

- [ [ ( [ Y^^=^+yYY^/j==dz\ (t>e{y)dyd: 
JtJt \Jt 2-sJ fN{z) J 

' f f ( f 

JtJt \Jt 




< 


1 

< - 
- 4 


1 

4 


CUN In 


{t, z)dz / ty^^x+y\{w)aN{t, w)dw 4>e{y)dydx 


I ^ \ \ f f f 

-{t,z)dz] / / / onY, w)li^_y^^]{x)(f>eiy)dxdydw 

J JtJtJt 

-{t,z)dz] ( / aN(t,w)dw) ( / y4>eiy)dy 
OiNjN / V4T / V-'T y 


OiNfN 

|V/iv|' 


By (14.221) . (14.231) and (|4.25p we obtain 


r f I/tv - (/tv)J < JcMTf [ yUy)dy 
Jo Jt y Jt 


(4.24) 


(4.25) 


(4.26) 


which completes the proof of the first part. 
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For the second part, the property ajv = 1 enable us to enhance the calculations of (|4.26l) in a 
way that 


l/v - ifN)e\^ dxdt 


Jt 

< [ [ x) - fN{t, X + y)f (j)eiy)dydx 

Jt Jt 

|2 


< 


4 /X 


^ (t>e{y)dydx 


4 / / / \z-y,z]{x)(t)e{y)dxdydw 

4 4T 4t JVC, Z) 

f f \VfN{t, 


/ 


4 Jt /v(t, 2 ) 




/ yMy)dy 

Jt 


and we are done. 


□ 


The following lemma is a summary of elementary convergence results which are useful in our 
context. 

Lemma 4.11. Let {fN}^=i, {9 n}n=i sequences of functions on [0, T] x T. 

(1) If In f, 9 n ^ 9 strongly in Li and II/atH^^ < C for all N, then f^g^ fg strongly in 

Li. 

(2) If fN ^ f strongly in L 2 and g weakly in L 2 then fN9N f9 weakly in Li. 

(3) Assuming that /at, / are weakly differentiable and g]\f, g > 0 for all N. If fN f weakly in 
Li) 9 n g strongly in L 2 and fN/9 n}'m=i uniformly bounded in L 2 , then ^fn/9N 

V f /g weakly in L 2 . 

Proof. (1) For any M > 0, 

f f l/vS-v - fg\ dxdt < f f I/atI \gN - 9\ + Is-I l/v - /| dxdt 

Jo Jt Jo Jt 

< f f C\gN - g\+M\fN - f\ + 2C\g\ligiyMdxdt 
Jo Jt 

then we can send N ^ 00 and then M ^ 00 to obtain the desired result. 

(2) For any bounded function U, 

f f (UfNgN - Ufg) dxdt 
Jo Jt 

< f [ \U\\gN\\fN - f\dxdt + [ [ {U fgN-U fg)dxdt 

Jo Jt Jo Jt 
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and since {gN}^=i is uniformly bounded in L 2 the first term converges to 0. Uf is a L 2 function 
and therefore the second term goes to 0 as well. 


(3) For any subsequence of {V/Ar/g'jv}?/=i) we can take a further subsequence which converges 
weakly in L 2 to some u. Then it suffices to show u = ^ almost surely. To this end, without loss of 
generality, we assume ^ u weakly in L 2 instead of its subsequence. Then, V^ QU weakly 
in Li by (2). However, for any smooth function v, 

/ / v{gu)dxdt = lim / / vVfj^dxdt = lim — / fj^Vvdxdt 

Jo Jt Jt N^oo Jq Jy 


f f fVvdxdt = f f vVfdxdt 

Jo Jt Jo Jt 


and therefore we obtain V/ = gu. 


□ 


Now we are ready to prove the lower semicontinuity of the rate function. 


Theorem 4.12. The functional is lower semicontinuous. 

Proof. It suffices to show that if x)dx p{t, x)dx weakly in (^([O, T], ./#(¥)) and ^ 

M for all k then I^i^{p) < M. Since we already assumed the initial LDP as in Assumption 2, it 
is enough to consider the dynamic part. 

We start by considering a functional 


x)dx - [ (j)^p{0, x)d: 

Jt Jt 




dxdt 


(4.27) 


on then 


^Tynip) = sup A^{p) 
(/>eC°°([0,T]xT) 


and hence it is enough to show lim^^oo = A^(p). The convergences of the first three terms 

in (I4.27P are direct from the weak convergence of p^^^ and therefore it suffices to show 

^ D{p)Vp weakly in Li([0, T] x T) (4.28) 

A{p^^^) —)• A{p) strongly in Li([0, T] x T). (4.29) 

First note that the uniform boundedness of I^ioriP^^"^) nnd Lemma 14.91 together imply 


// 




T P 


T {X + p{k)^p, 


(k) 

'c 


< M', Vc 


(4.30) 
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for some M'. Thus, satisfies the conditions of Lemma 14.101 with = A + and 

therefore convergence of pc is strong in Li. Moreover, satisfies the conditions of 

the second part of Lemma 14.101 because of Lemma 14.41 and hence p strongly in L 2 ■ 

(k) (k) 

To show (14.281). first note that the cth element of is ^ and 

therefore it suffices to show that 


Vpi^\ Vpc 

A+p(^) A+p 

- ^Vp 


weakly in Li([0, T] x T) 
weakly in Li([0, T] x T) 


(4.31) 

(4.32) 


for each c. Note that (I4.31|) follows directly from (14.301) and (3) of Lemma fd. Ill For (|4.32|) . by the 


same argument as before, we can show 


Vp(fe) 


Vp 


weakly in L 2 and it is also easy to check 


that 




Pc 


y/ X+p(k) y/A+P 

To prove (I4.29p . we need to show 


yT+pW vT+p 
strongly in L 2 - Thus, by (2) of Lemma [4. Ill we can prove (I4.32p . 


A±p« (fc) 
A+p('=) 

p''c'> Ik) 
\+p^k) Pc' 




strongly in Li ([0, T] x T) 
strongly in Li ([0, T] x T) 


(4.33) 

(4.34) 


\_L \_L 

for each c, c'. Since ^ ^ and 

prove (I4.33p . (I4.34p by (1) of Lemma [4.111 


strongly in Li and bounded by 1, we can 

□ 


4.3. Upper Bound. In this section, we establish the LDP upper bound for with the 

rate function The upper bound is usually based on the exponential martingale with a mean 

of 1 and the martingale should be suitably chosen such that it can be approximated by the density 
fields An exponential martingale such as this can be built by first using (t) (c/. (|3.15|) f in 

a way such that 

V E E 2 ,"(*))<iT(*). (4.35) 

where g = ••• , ™'^)[ G T] X T)™'. This martingale can be reinterpreted as 

r-T 


1 

N 


EE 

C=1L 


zf{T)) - g(^\0, zfm - [ g[^\t, zf^{t))dt 

40 


(4.36) 


according to Ito’s formula. As we commented in Section [321 represent zi^{t) — zj^{0) as 


A^A + 1 


^ V(aTI) E A(‘) + V(aTI) E («£(*) - M. 

^ ^ ^ k\k^i h-h-Ao 


(t)) 


(4.37) 


k'.k^i 
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and therefore the quadratic variation is 


A2 2A + 1 


dt -\- 


dAf{t) 


(A + 1) 


2-- -j 


(4.38) 


(A+ 1)2 iV(A + l)2 

Note that, although this expression relates to the local time, we can replace it by the local density 
by using Theorem 12.11 However, even after that, we still have a problem in (|4.36p . Broadly stated, 
we have a nuisance term relating to pt in the final stage that should not have appeared. The strategy 
for eliminating this term is to add another martingale 

[ Jxit, xf{t))dl3i{t) (4.39) 


i=l 


to (I4.35P with a suitably chosen J G C^’2([0, T] x T), which also has an alternative representation 
TV p '~p / \ 

J(r,xf(T))-J(0,xf(0))-^ (ji + ^J..j(t,Tf(t))dt (4.40) 


according to Ito’s formula. 

We now start the proof of the upper bound by defining a martingale J) for g G (^^’^([O, T] x 

T”*) and J G T] x T) by 

MN{g, J) 

'Y' Y' ^ 

= [ + f xf{t))dpi{t) 

do c=u^tn Jo i=i 


c=i 

.T - r 


rl 

L 

^ c=iie/+ 




dpi 


+ E E 

by (I4.37P where 


l<Ci<C2<mi^lN j^jN ■ 




- m 

GN{t, x^{t)) = ^ E E 


c=iie/A 


z^it))=gtHt, z^{t))-gtHt, z^it)). 

The next object to be characterized is A]sf{g, J) which must satisfy 

EArexp|MAr(5, J) - A^ig, J)| = 1. 


(4.41) 
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We can find such an A]\f{g, J) by 

AN{g-, J) 


1 pT m 

“ c=iieiN 


A 


^ jSf'(*. (t)) + x«(«)) + Mt, ref (()) 


1 2 


dt 


+ XN 


mx+i) ] 


{dAf^it) + dA^{t)) 


l<Cl<C 2 <m 


1 

2 

+ 


rT ^ 


EE 

c=Ue/A' '■ 

A 


A 




n 2 


dt 


2iV(A + l)2 


E E r E' 


dAfc2{t) + OAr(l) 


l<ci<C2<m ieL^ 
^1 


2 

where C/(x) =e^ — x—1~ Note that the error term is OAr(l) because of (j2.5ip . 

The next step is to approximate Mjv(5, J) and A]\f{g, J) by a density field of /r^(-)- To carry 
this program out, we define a set ^N,e,5 C ^([O, T], T^) such that x(-) E ^N,e,S if and only if 


[ x{t))dt < 6 and f x{t))dt 

io ’ Jo ’ 


<5 for c = 1, 2, • • • , m 


(4.42) 


where 


1 ^ r 1 

x) Zi) ^ - Xi) - ((i+(xj - Xi) + (5+(xi - Xj)) 


i=l 


j-j¥=i 


E E I 

l<ci<C2<m V 

ci’*^ C2 




2e 


Xeixj - Xi) - {6'^{xj - Xi) + (5+(xi - Xj)) 


Recall from Theorem [Q that 


limsuplimsup — Pat — —oo. 


r>n N 


(4.43) 


for any <5 > 0. Now, we can approximate M]\[{g, J) and Ajq{g, J) by the density field for x^(-) € 
t^N,e,5- More precisely, (|4.36p and (14.401) imply 

M^ig, J) = N [ch,,g,j(/i^(.)) + 0(<5)] (4.44) 
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for {■) G 3§N,e,S where the functional on C([0, T], ^(T)"*) is dehned as 

^e,g,j(^-) = (^T, giT, F^{T, x))) - (tto, ff(0, F^(0, x))) 


I 


n, (gt + gxx ) {t, F^{t, x))) dt 


2(A + 1)2 


+ (ttt, J{T, x)) - (tto, J(0, x)) - f 

Jo 

for TT. = (tt.^ 7 r. 2 , • • • , 7 r™)t, m = Y 1 T =1 and 


'^ti ( Jt F ^Jxx ) {t: x) ^ dt 


F^{t, x) = x + inidy), v{y - x)) 

Note that we used (I4.42p to replace the local time by the local density vr^ * Similarly, 
obtain 


AN{g, J) = N 


^ergAt^^\-)) + 0{5) + 0[- 


for x^{■) G ^N,e,5 where 


A 2 

= 2 {X + i )2 si {t, F^it, x))) dt 

A A 

F ^ ^ ^ J {.Jxid: x) F Kj^^g(t)) Qx {t, F-„^(t, x))) 

J {J^{t, x) F Ki-g{t)f'^ dt 

rT 

/ At, At * ie) {x) Ait, F^{t, x))) dt 
Jo 

/ At, L^,g,eit, x)gxit, F^t, x))) dt 

Jo 


dt 


F 


F 


A 

2(A + 1)2 
A 

2(A + 1)2 Jq 


with 


Ait, x) = (gAit, xf, gx\t, x)^ • • • , AAit, xf 
Ffr^At) ~ ^ At: gxit, Fj^[t, x))) 

Ln,gAt, x) = At * ie) it, x) ■ Qx (t, F^t, x)) . 

We can combine (I4.4ip . (14.441) and (I4.46P so that 


— logE^ [exp {A^ - 4'e,g,j(A (•))] } • 


(4.45) 


we can 

(4.46) 


(4.47) 


(4.48) 


1 
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Now we are ready to establish the large deviation upper bound for compact sets by the standard 
method (e.g., Chapter 10 of [l7]). For any open set 0 C (^([0, T], ^(T)"*), 

lim sup log Qn [^] 

N^oo JV 

< max| lim sup log Pat & 0} H ^N,e, 5 \ , hmsup ^ log Pat [■^n,€, 5 ] 


1 N—>-oo ^ 

and by Chebyshev’s inequality with (14.481) . 

■logP at £ 0} n ^N,e,S 


A^^cxd 


N 


1 


< ^ logPv [exp {N - d',,g,j(/2^(-))] } • 


N 

- {^ 6 ,g,j(T.) - d',,g,j(if.)} 

7T.£i/ 


= “+ 0(h) + 0 


Hence, we obtain 


hmsup ^ log Qn [^] < inf , sup ^eAN,g,j{Ti'-) 

N^oo 


where ^e,S,N,g,j{^-) is dehned as 

maxi - (d>e,g,j(7r.) - d'^,g,j(7f.)) + 0{6), lim sup log Pat 

L N^oo JV 

Then, by the Minimax lemma (c/. Lemma 3.2 of Appendix 2 of |17)L we have 


A,e,(5] 


limsup^logQAf [=^] < sup inf 0,^^s,N,g,j{^-) 
N^oo iv i.ex 


(4.49) 


for all compact sets C (^([0, T], Notice that, by Lemma 14.71 

lim sup ^ log Qat = lim sup ^ log Qn n ^J^ior] 

V^oo iV N^oo iV 

and therefore we can replace sup^ in (I4.49P by sup^^. . Moreover, for p G 

we havi^ 

^e,~g,j{p) - ^e,~g,j{p) = ^g,j{p) “ ^g,j{p) + 0 ^( 1 ) 


7 


The precise form is ^e,g,jipi-, x)dx) and so on. 
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where ^g^j{p) and ^g^j{p) are derived from ^^^g^j{p) and ^^^g^j{p) respectively, by replacing p * te 
and p* Lf:hy p and p. Thus we can rewrite (I4.49j) as 


lim sup ^ log Qat [TT] 

N^oo JV 


< sup inf max 


- (^g,j(p) - ^g,j{p) + Oe(l) + 0('5)) 


limsup^logPjv [^N,e,&] 

N^oo JV 


Now letting e ^ 0 and then <5 —?■ 0 so that we obtain, 
lim sup — log Qat [ < — inf 

N^OO N 

Consequently, it suffices to prove the following lemma. 
Lemma 4.13. For each p G 


sup 

g,J 


{d>g,j(p) -Tg,j(p)} 


(4.50) 


sup {^g,j{p) - ^g,j{p)} > iTynip)- 

geCU2([o,T],T)’" 

JecU2([o,r],T) 


(4.51) 


Proof. We first assume that p G C'^’^([0,T] x T)™. In this case, 

Fp{t, x)=x + / u{y - x)p{t, y)dy 

A + i Jt 

satisfies VT), = (c/. Proposition 5 of [10] ) and hence invertible for each t. Let Gp{t, x) be its 

inverse and then derivatives of Gp are given by 


d^Gp{t, x) = . 

A + p{t, x) 

dxxGpit, x) = 


(4.52) 

(4.53) 

(4.54) 


{X + p{t, x)) 

dtGp{t, x) = - ^ —- / i>{y - x)pt{t, y)dy 

A + p(t, x) Jj 

For given / = (/«, /(^l, ... , /M)t e Ci’2([o,r] x T)-, we take corresponding g and J by 

5(L x) =f{t, Gp(t, x)) (4.55) 

Tu ^ f YlT=iPc{t^ y)fx\t, y) 

J{t, x)= \ , — 77 — 7 - u{y-x)dy. 4.56) 

Jt X + p[t,y) 


Under these choices, we will show that 

^3,j(p) - ^g,jip) = A/(/5) 

where Aj{p) is defined in (I4.27jl . 


(4.57) 
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We first compute ^g^j[p). The main trick is to rewrite dtGp{t, x) in (|4.54p as 


~ic) ~ 


Px{t, x) 


A + pit 

and then we obtain 

^gA'p) 

= [ Pp{T, x)dx - [ /t,5(0, x)d. 
Jj 

-/ 7 " 


2{X + p{t, x)) 


(4.58) 


2{X + 2{X + p)^"" 


p{t, x)dxdt 


+ 


+ 


fL 


A + p. 


1 


v{y -x)[pt- -pxx ] {t, y)dy 


flKt, x)dxdt 


f J ■ p{T, x)dx — f J ■ p{0, x)dx — f f 
Jt Jt Jo Jt 


1 


Jt ~i~ ^Jxx 


p{t, x)dxdt. 


With our choice of J, the third and fourth lines cancel each other by the integration by part. 
To simplify ^g^j{p), we start from an observation that J satisfies 

Jx{t, x) = -Kppg{t) + x) 

A p 

where Kp^g{t) is defined in (14.471) . This enable us to compute 

^gAp) 


(4.59) 


fl 


A2 


2(A + p)2^ 


m „ , / m 

S(T>) ETV 

c=l ' ' \c=l 


+ 


p 


2{X + pf 
X 




,c=l 
' m 


2(A + p)^^ 


m 

E(/. 


(c) 


PPc 


2n 


c=l 


2(A + p)2 

r [lvpA{p)Vfdxdt. 

Jo Jt^ 




dxdt 


Now, (|4.58p and (14.591) complete the proof of (|4.57p . 

For general p G approximate p by p^. For given / G we can take and 

as (I4.55P and (j4.56p which correspond to instead of p. Then, by the previous step, we have 


^g„J.(Pe) - ^ge.APe) = ^f{Pe)- 
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Note that - ^~g„jAp) and are o,(l) since 

|Gp,(t, Fp{t, x)) -x\ = \Gp^{t, Fp{t, x)) - Gp^{t, Fp^t, x))| 

< x) - FpXt, x)| 

< jWpeit, ■)- p{t, OIIli 

for each t because dxGp^ = ^ Moreover, we have Aj{p^) —)• Aj{p) as e —>■ 0 as in the 

proof of Theorem 14.121 and hence 

sup {^ 9 ,jip) - ^g,jip)} > ATp). 
gecn2([o,r],T)"‘ 

Jecn2([o,T],T) 

holds for each / G and p G ^^lor- Consequently, we can complete the proof by taking 

supremum over /, □ 

Heretofore, we have established the large deviation upper bound for compact sets with the rate 
function but we can easily improve this result to the rate function where the argu¬ 

ment may depend on the initial conhguration. Moreover, since we have the exponential tightness 
by Theorem 13.11 the upper bound also holds for closed sets. 

Theorem 4.14. Under Assumptions 2 and 3, satisfies the large deviation upper bound 

with the good rate function and scale N. More precisely, for every closed set ^ G C'([0, T], 

we have 

limsup^logQiv [^] < - inf I^iorip)- 

N^oo ty p{-,x)dx£'^ 


4.4. Lower Bound. 


4.4.1. Perturbed Process. In general, the large deviation lower bound for the interacting particle 
system can be derived by observing the limit behavior of a suitably perturbed system. Stating 
this succinctly, the lower bound can be obtained by computing the relative entropy of such a 
perturbed process with respect to the original process. Thus, we should start by carefully defining 
the perturbations. 

Basically, we perturb our system in two ways. First, if the color of particle x^(-) is c, then we 
add the drift bfit, xAf)) to this particle. Note that the drift function depends on the color of the 
particle. Second, we change the jump rate between different colors. In the original process, we 
have a jump process MA (t) along the local time AA(t), which is the Poisson jump process with a 
constant intensity \N. We will also change this jump rate to \N + 7ci,c2(^) xf{t)) if the color of 
particles xf and xAf) are Cj and Cj, respectively. Then our perturbations can be summarized 
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by the m-dimensional vector b and an m x m matrix F, where 

h{t, x) = x), b 2 (t, x), ■■■ , bm(t, x))^ 

r(t, x) = { 7 ci,c 2 (i, x) : 1 < Cl, C 2 < m} 
both of which should satisfy the following conditions: 

(1) b and F are smooth. 

(2) 6(f, x) = 0 and t{t, x) = 0 for t E [0, ry] for some ry > 0. 

(3) f is skew-symmetric : 7^,02 = - 7 c 2 ,ci and 7c,c = 0. 

Remark 4.15. The third condition is not artificial in that changing 7 ci,c 2 and 7 c 2 ,ci by the same 
amount or the presence of 7c,c does not affect the dynamic of thus, we can assume the 

skew-symmetry of F without loss of generality. 


Let be the set of all (6, F) which satishes all of these conditions. For each (6, F) E a 
canonical way to describe the perturbed process is the martingale formulation. Indeed, we can 
understand this process by the measure on (^([O, T], T'^) such that for any / E C{G]\[), 

=f{x^{t)) - f(x^(0)) - ^ ^ Af(x^(s))ds (4.60) 

m 

bc{t, xf{s))Vif{x^{s))ds 

c=UeRNJ0 




fix'" {s))dA^Js) 






l<ci, C2<m 






where 

fix) = Dijf{x) - {\N + 'yci,c 2 )ifijix) - fjiix)) 

is a martingale with respect to the original filtration. The martingale also can be represented 

as (USD. Remark here that the rigorous existence and uniqueness of this perturbed process are due 
to Girsanov’s Theorem. 

The perturbed process P^'" is not too far from the original process Pjv in the following sense. 


Lemma 4.16. For each (b, f) E .^0 1 < p < 2, 

\r\P- 


^ log Eat 


dF' 


N 


(iPiv 


<c+0|- 


where the constant C could possibly depend only on p, b, F. 
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Proof. By Girsanov’s Theorem, 


dF 


b,r 

N 


dFj\f 


= exp < 




E 


V, 


(ci,C2) 




c=i ie/S 


1<C1, C2<m,ie/^, je/^ 


where 


bc{t, xf {t))d/3i{t) - ^ bl{t, xf {t))dt 

yic„C,) [dM^(t)+AiVd4^(t)] 

-[ lcuc2{t^ xf {t))dAfj{t) 

Jo 


under P^v- Since (6, T) G we obtain 


Eat 


dF 


.6,r' 


(iP^v 


< 


gCA(p(p-i)+o(i))^^g^p {Cp(p - i)7VAi^(r)} 


(4.61) 


for some constant C only depending on b and T. The last expectation can be controlled by (I2.5ip 
and we are done. □ 


An important implication of this lemma is the following corollary. 

Corollary 4.17. Theorems \2.1\ and \3.1\ are still valid under F^f^ instead o/Ptv for any (b, f) G ^o- 

4.4.2. Limit Theory of Perturbed Process. Let be the probability measure on ^([O, T],.y#(T)™') 
induced by {■) under the process P^^. Then, G N is a tight sequence because of Corollary 

14.171 Now, we can characterize all limit points of this sequence as the solution of a certain quasi- 
linear PDE. 

Theorem 4.18. Suppose that (b, f) G and ^{x)dx weakly as N ^ oo. Then, the 

support of any weak limit of , N is concentrated on the set of pft, x) G which is the 

weak solution of 

| = iv.[o(«va-v. 

with initial condition 7 (x). 

Let be a weak limit of G N. Then is concentrated on due to Lemmas 

ITTI and 14.161 We start by studying the limit of the uncolored empirical density which is no 

more the solution of the heat equation. 

Lemma 4.19. Let p{-, x)dx be any weak limit point of Then, p = Pc satisfies 

^ = iAp-V{i.p) (4.63) 


Mp)( 


P-yp 


(4.62) 
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in a weak sense. 


Proof. By the Ito’s formula, 
N 


1 

N 


^f{T, xf(r))-l 5 ^/( 0 , xf(o)) 

i=l i=l 

rp ^ ^ rj-i ^ 

^ £ |/i + ^C(i)/x + (t, xf{t))dt + Jfj^ £ fxii, xf{t))dpi{t) 


where c{i) is the color of particle x^ {■). Then (|4.63p is straightforward since the last term is 
negligible. □ 


Proof of Theorem \4.18\ The main machinery is again zf{t) in (I3.15p . However, we should be careful 


since (t) is not a martingale under but instead satisfies 


dz^ {t) =d.y^l'^ {t) + 


N / 


^ ^ ^ c=ljel^ 


(4.64) 


A 1 

+ xf{t))dt + tyxTT E^'=o,c(i, xf {t))dAf^^{t) 


(A + 1)^ 

where cq is the color of the particle x^{t) and (t) is the martingale given by (I4.37p which was 
just zj^ft) under P^r. For given 

~g={gW,g(^),... T] X T™) 

we can apply Ito’s formula such that 


- m ^ m 

77 E E T(T)) - 77 E E ^-"(O)) = 01 + 02 + 63 + 04 


(4.65) 


c=i ieiN 


c=i ieiN 


where 


^ E E z^{t))dt 

C=1 i(ilN 

®2 = / 4 E E ^xKt^ zf {t))dJtl^{t) 

do c=l i(zlN 

©3 = /" ^ E E i 

C=1 i(ilN 

rT -I r Ntj.\\ 

04 = ^ E E 3x\d^ ^i{x^{t))dt + ^ 


c=l 
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where 


ai{x^{t)) = 


1 


A 


iV(A + l) 


Y hit, x^it))dt + Y^hit, xfit))dt. 


k=lj(.lN 

We first claim that 02 is negligible since the order of the quadratic variation is 0{1/N). In the 
formula (I4.38P for the Brownian part is easy to compute. For the quadratic variation of 

the Poisson part, we only need to check 


[A^iT)] = E' 


N 


1 

N2 






< c 


(4.66) 


for some C where E^'" denotes the expectation with respect to To prove (|4.66p . let us define 
RNix) = dixi — Xj) for X G where giz) = ^ C'(T). Then, by Tanaka’s formula 

(MD, 


RN{x^{T))-RNihhQ)) + 


N t 


TiN - 1 ) 






^ Z] / it)) a'{x fit) - xfit))dt 

: / I *^0 

d/3iit) 




1 N pN 

= 4^ f 




Y a'ixf it) - xf it)) 


and we can check (j4.66p by simply taking the expectation. 

Now we substitute d ^ ™ ®3 (|4.38p and then apply the replacement lemma for the 

perturbed process (Corollary I4.17p to mollify the local times in ©3 and ©4 by local densities. By 
doing so, we obtain 


lim sup lim sup 
e^O N—,-00 


jbT 
■ N 




N 


Y Y {-nT.Y 1 


C=1 i€lN t 


c=i ieiN 


[af^ + + ji(a^^(i))#) it, zfit))dt 


> S 


= 0 


ji(x^(t)) = afx^it)) + Y 


"fcjkit. 


X. 


N 


k=l 


A+ 1 


it)) {k)i 

—pjh 


N 


it))- 


where 
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Now we represent all the terms as a function of the density fields of fj-Ni') and then send N ^ oo 
along the subsequence of N along which ^ weakly. Then, 

rT 


limsupQ^’" 




P ■ 


^g,p{T) - 3g,p(0) - f Rg,p,e{t)dt 

Jo 


> 6 


= 0 


(4.67) 


where 


x)dx, x))) 

^g,pAt) ={p{t^ x)dx, ^gt + x) 

Here, kp,, = k^% • • • , kl^V is defined by 


P.C 


x) 

^ Fpjt, x)) 

A + 1 


' « m 

/ b{t, y)^p{t, y)dy + Xhc{t, x) + ^7c,fe(t, x)pk * ie{x) 

k=l 


for c = 1, 2, • • • , m. 

The final step is to substitute g(t, x) = f(t, Gp(t, x)) where Gp = F~^ is the function defined in 
Lemma [4.131 Of course, this is possible only for p is regular enough. However for general p G 
we can use Peit, x) = f{t, Gp^{t, x)) instead and then send e —>■ 0 at the final stage to obtain the 
desired result as in Lemma 14.131 We will not repeat this procedure here. 

For p G we can compute various derivatives of g in terms of those of / by using (|4.52p . 

(14.531) and (|4.54p . Furthermore, we can explicitly compute (I4.54p by using Lemma 14.191 in a way 
that 


I®""’ = - x + Mt,^) L ~ 


A + p{t 

—- [ v{y - x) I l-Ap{t, y) - V \b{t, y) ■ p{t, y) 
x) Jf [2 1 

b ■ p{t, x) - /tp b{t, y) ■ p{t, y)dy 


A + pit 

Pxit, X 


dy 


+ 


2 (A + p(f, x)) (A + p) 

where we integrated by part at the last equality. By letting e —?> 0 at (|4.67p . we obtain 


pit, x)dx : / f ■ piT, x)dx — / /-^(O, x)d: 
JT Jt 


:x - Hi - Ho = 0 


= 1 


































LARGE DEVIATION PRINCIPLE FOR INTERACTING BROWNIAN MOTIONS 


61 


where 


H 


1 — 


ft + 


H 


2 — 


i:i,( 

m „T n 


A 


2(A + p) 


fxx - ^-^TT^—^fx] {t, x) • p{t, x)dx 


^bc + 'Yfjk=l Pkbk 
\ + P 


2{\ + pY 

m 

'yc,kPk 


x)dx 


By performing the integration by part, we can rewrite H 2 as 


f I 


A{p){b-jtp 


dx. 


(4.68) 


Note that we used the skew-symmetry of T here. This completes the proof, since Hi and H 2 
correspond to • [D{p)'Vp] and —V • A{p) — yf/5^ respectively in (|4.62p . □ 


4.4.3. Uniqueness and Approximation Proeedure. The lower bound computation, based on the limit 
theory of the perturbed system presented in the previous subsection, also requires the uniqueness of 
PDE (14.621) . Let consist of p satisfying I^iorip) < 00 then C t>y Lemma ITTl If we can 

prove the uniqueness of (|4.62p for the class of then we can directly compute the lower bound. Of 
course, the uniqueness of a quasi-linear PDE such as (I4.62p whose diffusion coefficient is not elliptic 
is hard to achieve at the desired level of generality. Instead, we establish a somewhat narrower 
uniqueness result, which should entail an additional approximation theorem. Thus, Theorem 14.201 
gives the uniqueness result and Theorem 14.221 provides the corresponding approximation procedure. 
We remark here that our methodology in the current subsection originates from and is similar to 
the methodology described in Sections 5 and 6 of |24) : hence, some details, especially related to the 
approximation procedure, are common to all of theses sections and will be omitted. 

Let a subclass of be the collection of p which is smooth on (0, T] x T, solves (I4.62p for 
some (6, P) G and satishes 

min inf pAt-, x) > e for some e > 0 (4.69) 

l<c<m {t,x)e[r],T]xT 

where p comes from the second condition of Then, we can state the uniqueness theorem as 
following theorem. 


Theorem 4.20. Suppose that u G is a solution of (|.^.6'il[ ) for (6, P) G with initial condition 
7 o(x) which satisfies IfvG is another solution of the same equation with the same initial 

condition, then v = u. 


Since the diffusion matrix is not symmetric, the usual technique based on the propagation of the 
Sobolev norm of h — ft is not available here. Instead, we examine the relative entropy of v with 
respect to u which requires v G Loo([i?, T] x T™). This boundedness does not automatically follow 
from the membership of and therefore, we require an independent argument to demonstrate 
this. 
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Let V = Wc where v = (ui, V 2 , ■ ■ ■ , Vm)K then it is enough to show v G Loo([r/, T] x T). 

First note that v is the solution of the heat equation in [0, r]] and therefore v{r], •) is a bounded 
function. In [r], T], we can add each coordinates of (|4.62p to obtain the equation for v: 


vt 


1 

2 


Av — V 



1 


Av + V{bv) 


where 

1 ™ 

b= - "Y] bcVc G Loo{[rj, T] x T). 

c=l 

Therefore, we obtain v G Lao{[rj, T] x T) from the following lemma. 


(4.70) 


Lemma 4.21. Suppose that w is the weak solution of 

^ = ]^Aw + V{hw) (4.71) 

with the bounded non-negative initial condition wo{x). Ifb£ Loo([0, T] x T), then w G Loo([0, T] x 
T). 


Proof. We first extend the equation to M. More precisely, we periodically extend 6 to M and call it 
b and then, consider the equation 


dw I ^ 

-^-Au, + W 


bw 


(4.72) 


where the initial condition is Wo(x) for 0 < x < 1 and 0 otherwise. To analyze (I4.72p . let us consider 
the diffusion 

dXt = dWt - b{t, Xt)dt 

on M where Wt is standard Brownian motion under the Wiener measure P. Note that the existence 
and uniqueness of Xt are guaranteed by Girsanov’s Theorem. Then (14.721) is the forward equation 
for Xt and therefore w{-, •) can be represented as 


w 


(t, y)= P(0, x; t, y)woix)d3 

Jo 


where the p{0, x; t, y) is the transition kernel of Xf. To compute this kernel, we assume that the 
Brownian motion under P starts from x and then consider a probability measure Q on C'([0, T], M) 
defined by 


dQ 

dP 


= exp 


b{s, Xs)dWs -I- f b^{s, Xs)ds 

^ ./n 


SO that Xt is a Brownian motion starting from x under Q. Then, 
P{Xt G [y,y + dy])=¥p 
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and therefore the kernel can be written as 
p(0, x; t, y) 

= (?(0, x; t, y)E'^ exp 

where 


(4.73) 


g(0, x;t, y) = 

7^,y a: 


[ h{s, z^,^y)dzfy[ P{s,z^’y)di 

Jo ^ Jo 

1 -(y-x)2 


: exp ■ 


2t 

is the standard heat kernel and {Zg’^, ^s}s<t Brownian bridge connecting x at time 0 

and y at time t under Q. 

Our aim is to estimate the kernel p by using (I4.73p . Observe that Zg’^ satisfies dZg’y = ds + 

dWg where {kFs}s<i is a Brownian motion under Q. Therefore, we have 




exp 


|-^ b{s, z^’y)dz^’y - P{s, z-^^y)ds^ < a\ai (4.74) 


where 


Til = M'S exp I -2^ b{s, Z^^’y)dWg - b^{s, Zf’y)ds^ 

^2 = E^exp|—2 J b{s, Z'. 


Z : 

: 


t — s 

2 TTz-w* /!_ -r-Jz-sT/-* ^7^ "^y 


It is easy to see that Ai is bounded by exp{r||6 ||^}. For ^ 2 , note that has an alternative 
expression = £(T:£l±^ W where i® another Brownian motion and hence, 

we can bound A 2 as 


A 2 < e^l^-2'lE^exp<! C 


1 T_ 


t(t — s) 


ds 


(4.75) 


where C could possibly depend on b only. Now we have to estimate the expectation in (I4.75I) . By 
Jensen’s inequality. 


E<3 


exp 




W. 


= E*3 


exp 


t(t — s) 

1 


ds 


(4.76) 


z2CViy/t — s 


< 


0 2y/t\ft — S 


W. 


t(t — s) 


W 


t(t — s) 


ds 


ds 


< 2e 


2C^T 


By (j4.74l) . (14.751) and (I4.76p . we obtain an estimate for the kernel p as 

p(0, x; t, y) < x] t, y) 


where constant Ci, C 2 only depend on b, T. This kernel estimates implies the uniform boundedness 
of w. □ 
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Proof of Theorem \4.20 Since the equation is linear parabolic with smooth coefficients on [0, ry], the 
uniqueness is automatic at there. Thus, it suffices to establish the uniqueness on [rj, T] x T. By 
Lemma I4.21[ we know that not only u, v but also Uc, Vc for all c are uniformly bounded by some 
number M > 0 in this region. We can define the relative entropy at time t such a manner that 

Veit, x) 


p nv 

Hit) = / 

C=1 


x) log 


Ucit, X 


-dx 


then by the elementary property of the relative entropy, 

m 


“ Vucit, x)| dx>^Kit) 


(4.77) 


where Kit) = {'^cit, x) — udt, t)}^ dx. Note that 77(?y) = 0 and therefore we can compute 

Kit) as 


Hit) = f [ 

Jri Jt 


d. 


rj JT 



where 


u 


Vl 


Tj Jt 


V2 


'^Vcis, x)log 

_C=1 

t 


log- 

u 


Vejs, x) 

Ucis, X) 

t 


dxds 


(4.78) 


dtv — ( — ) dtudxds 


fog w = log —, log —, • • • , log 


t 


and - = 
u 


Vl V2 
Ui U 2 ' 


Vm \ j 
Vm . 


Ul U 2 Vm, 

Now, we replace dtu and dtv by the RHS of (|4.62p and then apply integration by part. At this 
point, the only object that we cannot control is Vu and therefore we should simplify the result to 
the following form: 


[ f-WAVv 

Jr, JT 


+ B\r + Cdxds. 


If we carry out such a computation, then the result is given by 


1 
2 


+ 


HL 


8 ( 7 ) 2 Vu + 8(7) 2 


U-6 - G - iu+x(^t)V7 


dxds 


(4.79) 


8(7)- 


U'6-G + -U-x(«)V7 


dxds 


where 


8(7) 

G 


xiv)Aiv)±xiu)Aiu) 


x(7)A(7)x(7) 

T [x(^)^(^)f'7 


xiu)Aiu)ru 


Now, we will ignore the first term in (I4.79p . For the second term, note first that each elements of 
U“ and G are bounded by \uc — Vc\ for some constant C. Moreover |x(7)V7| is uniformly 
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bounded by (|4.69p and where 1^ is m x m identity matrix. Thus (|4.79p 

is bounded by C K(s)ds for some constant C. Thus, the uniqueness follows from Grownall’s 
Lemma. □ 


Since our uniqueness theorem is not for the class of but instead for we need an additional 
approximation procedure. Since the rate function is lower semicontinuous by Theorem 14.121 it is 
enough to establish the following theorem. 


Theorem 4.22. For each p G we can find a sequence C p^^\0, x) = 

p(0, x) for all k, p^^i p weakly and 

limsup/^„(p(^)) < iTynip)- 

k^oo 


In general, this procedure is not difficult if the rate function is convex. Unfortunately, within 
the context of our work, the rate function is not convex and requires careful analysis. For this 
purpose, we adopted the general method suggested in |24) . where comprehensive details can be 
found. Therefore, we only outline the whole procedure here; and additionally highlight selected 
points that do not directly follow from their result, due to the difference between our model and 
the SSEP. 

Our strategy is to divide the approximation into three steps as Jq” C C C where 
the two intermediate classes and are explained now. The subclass consists of p G 
that satisfies ^ [D{p)Vp\ for t G [0, rj\ for some r/ > 0. The membership of additionally 

requires that for some a > 0, Pc{t, x) > ap{t, x) holds for all x G T, 1 < c < m and t > p' for some 

0 < p' < p. 

The first step is to approximate by and which is Theorem 6.2 of |24) . The strategy is to 
estimate p G by pi'^i G dehned by 


x) 


R{t, x) 

< R{2p — t, x) 
J{t-2p, x) 


for 0 < t < ry 
for p < t < 2p 
for 2p < t < T. 


where R is the solution of dfR = • 

in [24] . 


D{R)VR 


with initial condition p(0, x). We refer the proof 


The second step is to approximate by and this step corresponds to the Theorem 6.3 of 
|24| . For this step, we first select a smooth increasing function e : [0, T] —>■ R satisfying e = 0 on 
[0, pi] and e = 1 on [ry2, T] for some 0 < ryi < p 2 < p- Then we can approximate p G by 


x) 



Pcft, x) + 


eit)pc 

k 


p{t, x) G (f™. 


One can find a proof of this step in |24| as well but we present a little bit simpler one. 
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Let US define 


^dyn 





OR 1 



-V • 

dt 2 

D{R)VR 


dt 




and then it suffices to show 


limsup/^„ [r]2,T]) < I^yn iP‘, [V2, T]) 

k—^oo ^ ' 


(4.80) 


since we can choose rp to arbitrarily close number to 772- It is easy to see that the rate function 
^dyni '! [^2) T]) IS convex on the set 

{ m 

--Y, Rc{t, x) = p{t, x) V(L x) E [0, T] X T 

C=1 

Since p^^\t, x) = (l — p{t, x) + x) on f > 772 where 

p{t, x) = {pip{t, x), p 2 p{t, x), ■■■ , PmP(t, x))I E Dp 



we have 

T]) < (1 - ^) ITyn iP', [P2, T]) + (p; [ 772 , T]) (4.81) 

due to convexity. We can easily check that 


lTynip-AV2,T])= [ 

J r?o 


r 

1 


Pi - 2^p 


dt < 00 


-1:P 


and therefore (|4.80p directly follows from (I4.8ip . 

The last step is to approximate by In |24) . this step has been carried out by Theorem 
6.4, which consists of Lemmas 6.5, 6.6, 6.7 and 6.8. In particular. Lemmas 6.5, 6.7 and 6.8 are quite 
robust and we can apply their arguments directly to our model as well. It would therefore suffice 
to show that a similar to Lemma 6.6 of |24) is valid for our model. This is verified by the following 
lemma. 


Lemma 4.23. Suppose that r and p are non-negative weakly differentiable funetions on T satisfying 


I 

Jt 


iVpp 


■dx < oo and 


|Vrp 


T P 


(A + p)r 


dx < CO 


(4.82) 


and r < p. Then, | | is a uniformly integrable family on T. 
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Proof. Notice that < 2 (^i + ^2) where 

/ vMy)iy 

Jt x/X + p[x + y) 


Te 

1 


A'?, = — 


n 2 


We can bound Ai as 

= - ( V 


/ . =\r{x + y)V4>eiy)dy 

JT\^/X + Pe{x) ^/X + p{x + y)J 


1 / \ 2 
1 /_ r 


•v/X-j-" 


< 


^ Tv ^ 

^ V VXTp 


< 


and hence this part is uniformly integrable by (I4.82p . 

By applying Cauchy-Schwarz’s inequality to A 2 , we obtain 


|V^l' I |Vp| 


\X + p)r 2Xp 


A 2 < 


x/X + pe{x) x/X + p{x + y) ^ 


Therefore, A 2 < j{Bi + B 2 ) where 

Bi = 

Bo = 


- A - Vp{x + y)) 


{^Uy)y 


e{y) 


, , ,{vuy))\, 

px + y)— , , . dy 

My) 

■dy. 



-ii/ 

i ((vT^))_ - Vp{^ + v)y 




Since fj dy = ^ for some constant C, we can bound Bi and B 2 as 

- J j (^x/P{x + z) - ■\/p{x + w)^ M^)M'>x)dzdw (4.83) 

B 2 < J^(^x/ p{x + z) - x/ p{x + y)j Mz)^^-^-^^dydz (4.84) 

respectively. Since G //^(T) by (I4.82p . we can conclude that RHSs of (14.831) and (I4.84p are 
uniformly integrable by Lemma 6.5 of |24| □ 

4.4.4. Proof of Lower Bound. Now we are ready to establish the large deviation lower bound for 




Theorem 4.24. Under Assumptions 2 and 3, {QA}^_r satisfies the large deviation lower bound 
with the rate function MiM')- other words, for any p G and its neighborhood 0, we have 


- MoIoAp) < liminf — logQiv 

N^oo iV 


(4.85) 
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Proof. Thanks to Theorem 14.221 it suffices to prove (|4.85p for p E For such a p, we can find a 
smooth function U{t, x) on [0, T] x T satisfying 

= iv . lD{p)Vp] - V [A{p)WU\ 

“■i 4”„(P) = 5 lo fj VI7tA(fi)V[f. 

We first assume that p(0, x) = po(x) so that dehne T,p^u by 

tp = VC/| . 

Then, by Theorems 14.181 and 14.201 we have 


= (6, f) E 


: b- 


lim eff) = l. 


N^oo 


for each (6, f) e ^p,u- Then we can estimate ^ log Qn [^] such that 


lim inf ^ log Qn [ff] > — _ inf lim sup 
N^oo N (fe,f)eSp,£/ N^oo 


N ^ dFN 


(4.86) 


(4.87) 


by the standard argument, e.g.. Chapter 10.5 of m- 

Now, we compute the RHS of (I4.87p . The first step is to recall Girsanov’s formula (I4.6ip to 
deduce 


1 


dP' 


, 6 ,r 

N 


N dPiv N 


y v. 

N Z-^ Z-^ ^ ^ TV 

c=l l<ci,C2<m 


(ci,C2) 


(4.88) 


where 


= [ bc{t, xf (t)) [dxf (t) - dif (t)j - ^ f blit, xfit))dt 
Jo ^ ^ Jo 


K- 


(ci,C2) 




L 




XN 




it) - f 

Jo 


7ci,c2(C xf{t))dyAt) 


and Jjj (t) is the jump process related with the martingale MA (t). In particular, under P^j 

dxf (t) = djSiit) + dAf it) + bcAt, xf it)) 

dJij it) = dMfj it) + (AiV + 7ci,c2(i, xfit))) dAf^ (t) 
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where ci and C2 are colors of particles xf {■) and Xj{-), respectively. Therefore, (|4.88l) can be 
rewritten as 


C=1 

+ ^ Y. [ 7ci,c2(*> + 0 


(4.89) 


ci<C2,iel'^ 

: defi] 

< 6 where 


To use the replacement lemma, we define a set 6) C C([0, T], T^) such that x(-) E 6) 

if and only if ^(t, x{t))dt 

x) 

~ Y1 7ci,C2(^’ Xi) 


l<ci, C 2 <m 


Y^Xe{Xj - Xi) - {5'^{Xj - Xi) + h+(Xi - Xj)) 


Then, by Corollary 14.171 ^n{^, is super-exponentially negligible and hence 


lim sup 


N^OO 


1 


log 


dF 


b,r' 

N 


N dFj\j^ 


= lim sup lim sup 
e^O N—>‘00 


1 


<5) AT 


dP' 


N 


N dFj\f 


On h), we can approximate (I4.89P by 


1 


+ 


YY [ xf{t))dt 

T / 

Y [ 7ci,c2(i> xf {t))p\'^^\xf {t))dt + 0 { 6 ) + o(^ 


c=l is/' 


1 


ci<C2,ie/jY 

Consequently, we can conclude from (I4.90D and (I4.9ip that 


lim sup E^^ 


N^OO 


log 


jbT 

■ N 


N dPjv 




ll,c 2 pciPc 2 {t, x) } dxdt 


since p{t, x)dx by Theorem 14.181 and 14.201 


(4.90) 


(4.91) 


(4.92) 
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To complete the calculation of the RHS of (I4.87p . we optimize (|4.92p over {b,T) € '^p,u- This 
can be done by the Lagrange multiplier method and the optimizer turns out to be 


- A 1 J 

br = ^-V?7c + ■;—;—■ ^ ^ Pk^Uk 


7ci,C2 — 


A + p 
A 

A + p 


^ + Pk=i 
- VC/eJ. 


With these optimizers, the RHS of (|4.92p becomes ^ fjVU^A(p)V[/ = 

By following this approach, we completed the proof when p(0, x) = p^{x). The case for the general 
initial condition is also easy to obtain by the same argument by tilting the initial configuration 
appropriately. □ 


We conclude this section by summarizing the results that were obtained for the LDP for the 
empirical density of colors. 

Theorem 4.25. Under Assumptions 2 and 3, satisfies the LDP with the good rate func¬ 
tion and scale N. In other words, for any measurable set A C ^^([0, T], we 

have 

- >1 iZiori^-) < liminf ^ logQAr(H) < limsup ^ logQAr(H) < - inf 
neA° N^oo N N^oo SeA 


5. Empirical Process 

5.1. Propagation of Chaos. We start by explaining the relationship between the propagation of 
chaos, which is the LLN of the empirical process, and the LLN of the empirical density of colors in 
a more general set up. 

Consider the empirical process Rn = ^ induces a probability measure P/v on 

.^i(C'([0, T], T)). The limit theory for {P/v}^_^ can be obtained by verifying the tightness and 
identifying the unique limit point. The tightness can be demonstrated by the general technique 
introduced in |26| . Regarding the identification of the limit point, the limit theory of the empirical 
density of colors plays a significant role. Suppose that the limiting particle density p{t, x) is the 
unique solution of a certain parabolic equation dtp = ££p with the initial condition /0®(x) under 
Assumption 1. Furthermore, assume that if we color the particles by an arbitrary number of colors 
such that Assumption 2 holds, then the limiting particle density of each color c denoted by pc 
evolves as the unique solution of the parabolic PDE dtPc = -^p Pc with the initial condition p^fidx), 
where is a time-inhomogeneous generator that could possibly depend on p{t, x). 


Remark 5.1. For our model, and s^p is given by (jl.l2j) . 
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Under these assumptions, we can compute the limit of finite dimensional marginal densities of 
the empirical process. For instance, we can calculate the limiting joint density 

N 


lim Ejv 
N^oo 


|{i : xf(0) € A, xfjt) £ B}\ 


N 


= lim Etv 
N^oo 


^lyl(xf(0))lB(xf (t)) 


_ 2=1 


(5.1) 


in the following manner: we color the particle xf{-) by color 1 if x)^(O) G A and by color 2 otherwise. 
If /r'^(O) ^ p^{dx), then p^{0) l.A{x)p^{dx) and therefore, we can compute the limiting particle 
density pi{t, •) of color 1 at time t by the solution of dtpi = Pi with the initial condition 
\A{x)p^{dx). Therefore (15.11) can be computed as /Oi(t, x)dx. We can use the same method to 
compute the joint distribution for any finite number of times, (see [261132) for details.) 

Therefore, any limit points of should be the diffusion process with the generator £/p. 

Consequently, we can establish the limit theory of {Pn}]^=i as soon as the uniqueness and existence 
of such a diffusion process with starting measure p^{dx) are valid. This general theory can be applied 
to our model if the initial limiting particle density is bounded. 


Theorem 5.2. Suppose that (0) p^{x)dx weakly for a bounded function p^{x) on T and let 
p{t, x) be the solution of the heat equation with initial condition p^{x). Then P]\f ^ 6p weakly where 
P is the unique diffusion process on .^i(C([0, T], T)) with the time-inhomogeneous generator s/p 
defined by M.PA) . 

Proof. The tightness of is a consequence of Theorem 13.11 and the limit theory for the 

empirical density of colors is presented by Theorem 14.11 The uniqueness result for the diffusion 
with the generator £/p and the starting density p^(x) which is bounded can be found in Theorem 4 
of do]. □ 

Even though we have suggested a proof of Theorem 15.21 by using the empirical density of colors 
as an intermediate tool, this result was already established in |10| in a different way. The stronger 
result in uni showed the diffusive scaling limit of one tagged particle to be the diffusion with the 
generator £/p and also showed that any two tagged particles are asymptotically independent. Of 
course, these results imply the propagation of chaos for our model. 


Remark 5.3. For the general starting measure p^{dx), our methodology is still valid for the tightness 
and the identification of the limit point step. However, the uniqueness of the diffusion process with 
the generator £/p causes a problem. If p^{dx) is a singular measure, then the uniqueness generally 
does not hold. However, HU suggested a way to circumvent this pathological phenomenon by, 
roughly speaking, appropriately decomposing each mass at a point into a left and right mass. We 
were also able to extend our result to this regime. 


The remaining part of this article is devoted to explaining the LDP corresponding to Theorem 
15.21 under Assumption 3. A methodology for the SSEP for d >2 has been developed in |24||^ and 

®The original result was valid only for d > 3 but extended to d = 2 in [20] . 
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relies on the LDP for the empirical density of colors and Dawson-Gartner’s projective limit theory. 
The robustness of their method is such that we can almost apply it directly to our model. The only 
thing that has to be checked for our model is a certain class of martingale problems 


5.2. Martingale Problem. When we define the rate function J’ (ff) for the LDP of the empirical 
process in the next subsection, what we need is the perturbed diffusions with the generator + 6V 
for an appropriate class of b. The existence and uniqueness of such diffusions are not trivial and 
should be proven independently. In this subsection, we carry this out with the help of the results 
in [25]. 


Suppose that p{t, x) is weakly continuous in time, weakly differentiable in space and also satisfies 

- (Vp)2 


/ p(0, x)logp(0, x)d: 

Jt 


X < oo and 


II 


-dxdt < oo. 


'0 Jt P 

Then we define a class consisting of measurable functions b{t, x) on [0, T] x T such that 


^ = -A/9 — Vibp) and 
dt 2 ^ ^ 


fl 


opdxdt < oo 


(5.2) 


(5.3) 


where the first equation is weak sense. 

For measurable function c{t, x) on [0, T] x T, we define the generator £/p^c by £/p^c = + cV so 

that 


^p,c 


A / (2A + /9)V/9 

^(AT^'^+V" 2(A + p)2 


V. 


(5.4) 


Then (|5.3I) implies that p satisfies dtp = ^p^p for each b G ^p. Basically, we want to build a 
unique diffusion process with generator £/p^b for b G ^p with marginal density p to define the rate 
function of empirical process. However, the coefficients of the generator s^p^b only have limited 
regularities and therefore the existence and uniqueness in the spirit of Stroock and Varadhan is not 
valid here. Although there are some results on general coefficients (e.g., HD), these usually assume 
uniform ellipticity for the generator. In our case, the diffusion coefficient is 2 {\+p) ’ may not 

be uniformly elliptic since p can be unbounded in general. For the SSEP, Quastel and Varadhan 
|25| solved this difficulty by limiting the sense of the martingale problem in a suitable fashion. They 
considered the solution of the martingale problem not to start from a specific point x but from some 
initial distribution Pq{x). By doing so, they achieved a proper existence and uniqueness result in 
this context. Of course, we shall follow their approach and the main result can be stated as follows. 


Theorem 5.4. Suppose that p satisfies \5.2J and SSp (j). 

(1) For each b G Mp, there exists the unique diffusion process on T with the generator s^p^b 
with the marginal density p{t, x) at each time t G [0, T]. 

(2) For eaeh measurable function on T satisfying 0 < R^f) < 0/3(0, •), there exist unique 
diffusion P^q with the generator s^p^b o,nd the marginal density R{t, x) which is the unique 
solution of dtR = ^pb^ with initial condition Rfi and satisfies 0 < R < Cp on [0, T] x T. 
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The proofs are identical to those in Section 5 of [25]. The only obstacle when we apply the 
argument of |25) is the fact that our model possibly has an unbounded density p{t, x) whereas the 
SSEP has an a priori bound 1. We can overcome this by proving the following lemma as a substitute 
to Theorem 3.12 in |25| . 


Lemma 5.5. Suppose that p satisfies 115.S\} and fi. For eaeh measurable function on T 

satisfying 0 < R^{-) < Cp{0, •) for some constant C, there exists the unique weak solution R{t, x) 
of the forward equation 

f) /? 

^ = KhR ( 5 - 5 ) 

with initial condition R^ and satisfying 0 < R < Cp on [0, T] x T for some constant C. Moreover, 
R also satisfies the energy estimate 


fl 


(Vi?)' 


-dxdt < Cl + C2 


fl 


iVpp 


dxdt 


0 Jt P 


(5.6) 


'0 Jt (^ + p)i? 
for some constants Ci, 02- 

Proof. For each b E Sdp, we define R = Then, it is easy to see that R E F^p,^ and is the 

unique weak solution of dtPe = 32^* h^Pe- Note that we can write explicitly as 




u = V 


A 


2(A + Pf) 


Vtt + 




2(A + Pe 


-R\u 


(5.7) 


and it is easy to check that this generator satisfies the conditions of Theorem 3.12 in |25| . namely, 


fl 


(Vp.)' 


A 


0 Jt Pe 


2(A + Pe) 


dxdt < 00 


fl. 


Vpe 


2(A + Pe) 


-R 


2(A + pe) 


A 


Pedxdt < 00 


since we have 


{^Pe? ^ gyp) 
~\p 


Pe 

i.e\2 


rn'^Pe = ^ < {b^pfi 


(5.8) 


(5.9) 


and Pe is uniformly bounded by some constant M^. Therefore, we can apply Theorem 3.12 of [25 
such that there exists a unique solution i?^ of 

5i?^ 


(5-10) 

with the initial condition i?e(a)) that satisfies 0 < i?*^ < Cp^ on [0, T] x T as well as the energy 
estimate 


fl 


{VR^ 


T (A + Pe)i?' 


< Cl +C2 


fl 


(VpY 


dxdt 


T P 


(5.11) 


for some constant Ci, C 2 . This energy estimate can be derived from (3.26) of |25| and (15.8p . 
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Our aim is to send e to 0 in (15.101) in a proper way. To this end, let us first prove that 






e>0 


and 


are uniformly bounded in T2([0, T] x T), respectively. The boundedness of the first of these terms 
is obvious and that of the second term follows directly from (15. lip . For the last term, since p £ L 2 
by dSJl), ^ ^ ^ 

f f dxdt < f f p^dxdt < f f p^dxdt 

io JT io JT JO JT 

by Lemma [4.41 Therefore, we can take a subsequence {ek}^=i which converges to 0 and also satisfies 


R^f^ R, 




A + Pek 

weakly in L 2 for some R, U and V, respectively. 
We now claim that 


[/ and 


V P<^k 


V 


U= and V =—. 

X + p y/p 


For U, we know that A + —)• A + p strong in L 2 by Lemma (4.101 and — is uniformly bounded 
in L 2 by (15.lip . Therefore by (3) of Lemma [4.111 we can verify that U = For V, by (2) of 

Lemma 14.111 we have ■ y/P<^k V ■ yfp weakly in Li and therefore V^fp = i? or equivalently 

V 

F= 4. 

y/p 

These weak convergences in L 2 imply that 

Vi?"'' 




_^ Vi? 

A + Pet, A + p 
Vp 


-i?"'' 


A + Pet 

Jj^kR^k 


X + P 
bR 


R 


(5.12) 

(5.13) 

(5.14) 


weakly in Li also by Lemma 14.111 More precisely, (15.121) is derived directly from our definition 

of and (15.131) holds because - > strongly in L 2 , due to the uniform integrability 

of the form of (|5.8p . Similarly, (I5.14p is obtained as a consequence of (2) of Lemma 14.111 since we 
have 6"'=.y/p^ —)• b^/]) strongly in L 2 by (15.9p and ■^= ^ weakly in L 2 as we observed before. 
Now, (j5.12p . (|5.13p and (15.141) allow us to take the limit in (|5.10p along the sequence and 

by doing so we obtain ^ = £/*^R. Consequently, we proved the existence. 

The energy estimate (15.6p can be obtained by repeating the argument of Theorem 4.1 in [25] . 
Although this theorem requires the L^o boundedness of R, our bound R G ^2(0, T, Loo(T)) turns 
out to be sufficient for applying their argument to our specihc diffusion coefficient 
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Finally, let us consider the uniqueness issue. Suppose that u, v are two solutions then we consider 
the evolution of which is a well-defined function since 0 < u, v < Cp, such a manner that 


/ - {s, x)dx — / -(0, x)(ix 

Jt P Jt P 


ff 

Jo Jt 


(u — V) / 1 

^-’vi-Vp-bp 


u — V „ 

+ 2 -V 


A 


P L2(A-|-p) 


V{u — v) + 


Vp 


2(A -I- p) 


— b \ {u — v) 


dxdt 


fl 


V 


(u — vY 


-2V 


fl 


u — V 


Vp — bp] dxdt 


A 


P J + p) 


V{u — v) + 


Vp 


2(A -I- p) 


— b] {u — v) 


dxdt 


A 


'0 Jt P^(^ + p) 

This computation guarantees the uniqueness. 


[(tt — u) Vp — pV {u — v)f dxdt. 


□ 


5.3. Large Deviation Theory of Empirical Process. We start by defining the rate function 
for empirical process. 

Definition 5.6 (Rate function for empirical process). Let Q E .y#i(C'([0, T], T)) has the marginal 
density q{t, x) which satisfies (j5.2p . JJSq Y Y FT [Q |P^] < oo for some b E where the diffusion 
is the one defined in Theorem 15.41 Then, we can fincj^ bg E tJ^q such that the corresponding 
diffusion process with marginal density q{t, x) satisfies 


pQ 


1 6 r 

/ (p{t, x{t))dx{t) = / (p{t, x{t))dx{t) 

Jo J Uo 


for any smooth (p (cf. Theorem 7.3 of |24)). Then, the dynamic rate function ^dyn{Q) is defined by 


^dyn[Q) = H 


Q 


P^Q 


+ 


\ ! f bgQdxdt. 
^ Jo Jt 


(5.15) 


For all the other cases, .^dyniQ) is defined to be infinite. In addition, due to Assumption 3, we 
define the Sanov-type initial rate function J^initiQ) by 

^(0, x) 


AnitiQ) = / ?(0, x) log 

Jt 

Finally, the full rate function is defined by 


p0(x) 


-dx. 


J{Q) = JdyJQ) + ^iniJQ)- 


®Alternative way to define bg is the unique minimizer of the relative entropy H [Q |T^] over b G PSq. 
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The functional J^(-) defined in this manner is lower semicontinuous and has compact level sets(c/. 
Theorem 7.4 of |24|). 

Now, we can state the LDP for the empirical process in a concrete form. The following theorem 
can be proven by the general method presented in Sections 7 and 9 of [24| . which relies on the LDP 
for colored system and Dawson-Gartner’s projective limit theorem (c/. Theorem 4.6.1 of [3]). 

Theorem 5.7. Under Assumption 3, {Pn}n=i satisfies the LDP with the good rate function J^{-) 
defined in Definition \5.b\ and scale N. In other words, for each measurable set A C ./#i(C'([0, T], T)), 

— inf ^(Q) < liminf 4-log Pat (j 4) < limsup 4-logPAr(^) < — 

Q&A° N^oo N N-^oo N Q£A 

where the topology is the usual topology of weak convergence for measures. 
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